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Volume reflecting materials constitute a class of dielectrics which are capable 
of o high level of internal scattering such that a significant fraction of an Incident flux 
emerges from the Incident flux surface. As such the primary modes of energy transfer in 
these materials are conductive and radiative transfer. The research performed under this 
grant has focussed on the calculation of radiative transfer in highly scattering materials, 
both by exact and approximate means^and the interaction of radiative transfer in such 
moterials with the conduction mechanism, in both steody and transient modes . 

An approximate analysis of radiative transfer in highly scattering materials was 
developed based on the Kube!ka-A4unk differential equations — a set of two differentia! 
equations representing the spatial rate of change of radiative half-fluxes within the 
scattering media. These opproximate solutions of the Kubelka-AAjnk equations together 
with analytic solutions for the steady state temperature distribution for two types of 
boundary condifionsare given in reference 1. These solutions show the influence of back 
surface reflectance, scattering power,incident rodiative flux parameter and boundary 
conductive flux parameter on overall reflectance and temperoture distributions. This 
radiation Field analysis, adapted to spherical geometry, was applied In reference 2 to the 
evaluation of the thermal performance of teflon and fritted quartz as heat protection 
^ materials for entry into the atmosphere of Jupiter, 

*The NASA Technical Officer for this grant is Dr. Phillip R. Ncchtsheim, 
NASA Ames Research Center, Moffett Field, California 94035. 



Two exact analytic solutions for transient temperature distributions in non-emitting, 
plane pcrallel, diffuse reflectors were developed based on the approximate radiation field 
model presented in reference 1. The first of the initial-boundary value problems considered 

I 

(reference 3) deals with transient development of the steady state solutions for the tempera- 
ture in highly scattering media presented in reference 1. The solution proceeds from an 
initially uniform temperature distribution with an instantaneous change in temperature at 
one boundary and a zero conductive flux at the other boundary. The theory is applied to 
the transient heating of a 1 cm. quartz slab in reference 3. Results obtained are similar to 
those presented in reference 2 for the more complex case of cylindrical geometry. However, 
the solution is much simpler for the plane parallel geometry. Reference 3 also presents an 
evaluation of the adequacy of the Kubelka-AAunk reflectance solution for zero absorption 
when opplled to the determination of the reflectance of highly scattering and weakly 
absorbing media. The importance of using accurate overall reflectance values In the 
solution is demonstrated . The influence of rear surface reflectance, scattering power, and 
absorption coefficient-scattering coefficient ratio on overall media reflectance Is also 
briefly discussed based on the Kubelka-^VJnk theory in reference 3. 

the second analytic solution referred to above treats the more involved un- 
steady temperature development arising from specified constant radiative and conduc- 
tive fluxes at one boundary and vanishing temperature gradient at the other boundary 
(reference 5). It was shown that the solutions increase or decrease in time monot- 
onlcaliy for a given incident radiative flux depending on the volue of the boundary 
conductive heat flux. A relation was also obtained defining a singular condition for 
which a steady state exists. This relation obtained by study of the analytic solution 

may also be obtained by means of a steady state energy balance on the diffuse 

\ - 



reflector. Sample results showing unsteady temperature development for cases of asymp- 
totically increasing and decreasing temperature and for the critical case are also 
presented in reference 4. Some of the details of the^analyses not presented in references 
3 and 4 are given in reference 5. . 

While approximate solutions such as those discussed above are useful in making 
rapid estimates in preliminary design and in guiding qualitative thought, more precise 
representations of radiative transfer and temperature fields may be required for design 
purposes. Toward this end a method of coupling solutions of the equation of radiative 
transfer and energy equation was developed (reference 6). This complex technique, 
utilizing the method of idempotenfs to generate a quasi-steady radiation field, determines 
an instantaneous temperature field by solving the thermal energy equation (including 
radiative flux divergence term) using an explicit finite difference scheme . The resulting 
temperature field then couples with the radiative transfer solution through radiative emis- 
sion. The method was applied to the analysis of the transient development of the coupled, 
steady state temperature and radiation field distributions between two opaque, partially 
reflecting boundaries. The Intervening media were allowed participation in the energy 
exchange through the mechanisms of radlotive absorption, emission and isotropic scattering; 
and through the mechanism of thermal conduction. Reference 6 shows that this method 
accurately reproduces transient, non-scattering results in the literature as special cases. 
The method was also used to generate steady state results for comparison with existing 
Isotropic scattering solutions in the literature. Reference 7, an extension and refinement 
of the work of reference 6, presents previously unavailable transient solutions Tor the 
plone-parallel radiation field problem In which isotropic scattering Is present. These 
parametric studies of the influence of the effects of optical thickness, albedo, boundary 



emissivity and conduction-radiaHon parameter on the temperature and energy flux 
distributions in semi-transparent media may be regarded as exact, within the limitations 
of the spatial and directional discretization Inherent in the explicit finite difference 
representation of the energy equation and the Gaussian Quadrature representation of 
the transfer equation scattering Integral. The very complex computer program 
developed to obtain these solutions was given the acronym. CURCES for the Combined 
Unsteady Radiative and Conductive Energy System. A listing of the program is given in 
an appendix of reference 6, 

While the CURCES program provided the milestone In unsteady coupled radiati ,** 
and conductive energy transfer published in reference 7, it was cumbersome to use and 
lacked flexibility for adaptation to more complex problems. For instance, the idempotent 
method used in CURCES is mathematically incompatible with spatially varying radiative 
property distributions. As a result further solutions of the radiative transfer equation were 
obtained through the develop ment of a series of programs using an iterative technique and 
therefore bearing the acronym ITERAD for Iteration of Radiation. A discussion of the basic 
solution method and its convergence is given in reference 8. A comparison of the reflect- 
ance of a specific diffuse reflector as computed by the CURCES and ITERAD programs is 
given in Table 1 , 


Number of 

TABLE 1 

Reflectance* 


Ordindtes 

CURCES 

ITERAD 

2 

.8678 

;8672 

4 

.8649 

-8646 

6 

,8637 

.8633 

8 

.8632 

.8627 

10 

-- 

.8625 

12 

— 

.8623 

14 

— ' 

.8622 


*Albedo = .9995, Optical Thickness 3.177, Rear Surface Reflectance = 0.8 


It IS clear that the two sets of results agree very well and demonstrate the improve- 
ment in accuracy and diminishing returns associated with increasing quadrature order. 

Tliese results may also be compared with the value .875 obtained from the Kubelka-Mjnk 
equation (5) of reference 1 . Further evidence of the consistency of the CURCES and 
ITERAD programs is shown in figure 2 of reference 8. There individual radiation field 
intensify distributions inside a highly scattering medium calculated by the two methods 
are shown to be in very good agreement. Reference 8 also shows good agreement of 
ITERAD reflectance calculations with values from the literature over a wide range of 
radiation parameters. 

Through adaptation of the ITERAD program a study of the influence of anisotropy 
on the reflectance and internal radiation field of a highly scattering material was made. 
Some results of this study for the phase function: P (9 J =w[ I +xcos 9] are given in refer- 
ence 8. Here uj Is the scattering albedo, 9 Is the angle between Incoming and outgoing 
beams and X is an anisotropy parameter such that one obtains ... 
net backward scattering for -1.0<x<0 

isotropic scattering for x = 0 

net forward scattering for 0 <x 1 .0 

The results demonstrate that net forward scattering allows the penetration of radiotion to 
greater depths than isotropic scattering and -causes ' greater internal energy conversion 
to thermal energy through media and rear surface absorption. The reverse effect Is obtained 
for net backward scattering. As might be expected backward scattering increases overall 
reflectance while forward scattering decreases reflectance. Reference 8 shows the 
influence of anisotropy of scattering on the radiative flux and radiative flux divergence 
distributions as well as on the intensity field. 



An extensive effort was put forward in this program to evaluate two flux theories 
with respect to the transfer equation . It Is demonstrated In reference 8 that the Schuster- 
Schwarzschild two-fiux equations may be obtained by integration of the equation of 
radiative transfer over two hemispheres for which the radiation field Is given by 



I V) j 0 
i"(T), ^ < 0 


It is shown that the Kubeika-Munk equations are identical to the Schuster-Schwarzschild 
equations when a simple set of relations exist between the scattering and absorption co- 
efficients of the two theories. It ts pointed out in reference 8 that the relations between 
the coefficients of the two theories depend on the nature of the intensity field, A case in 
the literature is cited for which a different set of relations is obtained. 

Computations are presented in reference 8 which compare transfer equation solu- 


tions with reflectance, radiative flux, and flux divergence obtained from Kubelka-Munk 
analyses. It is demonstrated that the Kubelka-Munk two-flux analysis is capable of 
reasonable approximation of transfer equation results. The approximate radiative flux 
and radiative flux divergence relations of reference 1 are shown to be accurate only for 
albedo very near unity. The Influence of the relation between the transfer equation 
parameters and the Kubelka-Munk coefficients on the radiative calculation comparison 
Is also examined. 

The iteration procedure for solving the transfer equation was applied to the 
determination of the radiation field and related parameters fora scattering medium with 
non-unity refractive Index, A detailed study (reference 9) showed that certain calcula- 
tions in the literature had significant errors because of the poor opproximation of the 
scattering integral with the standord Gaussian quadrature. There It was shown using an 
-approximate model and verified using transfer equation solutions that direct application 



of Gaussian quadrafure to the evaluation of the scattering integral can result in large 
errors in certain cases. This occurs because of the presence of a discontinuity at the 
critical angle for total infernal reflection . Application of the model also showed that 
the error for a given quodrature order is not necessarily a monofonic function of the 
index of refraction and that it is possible under certain conditions to increase rather 
than decrease error by increasing quadrature order. Reference 9 also shows the signi- 
ficant error reduction possible by tailoring the quadrature approximation of the scatter- 
ing integral to the critical angle for total infernal reflection . The model of reference 9 
is also briefly applied to demonstrate that no error exists in the scattering integral as a 
result of the discontinuity at [X =0 in the unity refractive index rudiation field. 

A study of several citematives for application of Gaussian cpjadratures to avoid,,: :., 
the discontinuity error pointed out in reference 9 is reported in reference 10. There, by 
extrapolation of dn Increasingly accurate sequence of transfer equation reflectance values, 
precise values are obtained which are employed to evaluate the adequacy for dtsconfinuity 
error avoidance of several (Russian quadrature combinations involving sixteen total 
directions, it was found that the combination of fifth order quadrature in each of the 
forward and rearward critical cones and sixth order quadrature outside yields the most 
accurate reflectance solutions for sixteen total directions (on the average). As a result 
the transfer equation iteration program was revised fo Incorporate this combination. A 
number of solutions have been obtained with the revised program for refractive indices 
of T.2 and 1.4. Reflectance values for layers of finite optical thickness and refractive 
indices of 1 .2 and 1 .4 ore tabuloted in reference 10 along with o comparison of unify 
refractive index reflectances from several sources. This collection is the first compari- 
son of several sources of unity refractive Index reflecfances>and the first tabulation of 



non-onify refracfive index reflectances for the conditions considered/which is known to 
the principal investigator. 

A final version of the tronsfer equation iteration program was developed which 
incorporated all the features in prior iteration programs as discussed above as well as 
certain other useful features. These feahjres include a capability of treating composite 
plane layers with several differing values of absorption and scattering coefficient and 
a capability of treating radiation problems with non-gray incident flux distributions and 
property distributions. This program called the Band Model Program is discussed in the 
accompanying appendix. A Usftng of the program and several test cases demonstrating 
the band model capability and the spatial property variation capability are included 
there , . 
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APPENDIX 

fINAL BAND MODEL PROGRAM DESCRIPTION 
By Dantel W. Drago 

Over the past fev/ years at the University of Tulsa, a program has been designed and 
developed to approximate the solution of the transfer equation as It applies to a scattering, 
absorbing, emitting, plane-parallel, semi-infinite slab subject to a diffuse, incident, radia- 
tive flux.* The slab is bounded on the rear by a specularly reflecting surface which is 
modeled using either uniform or Fresnel reflectance. The temperature distribution in the 
slab may be an arbitrarily specified function or constant . This program does not analyze 
changes in the temperature field over time but could easily be modified to do so. 

The incident radiative flux and subsequent solutions of the internal radiation field 
may be divided into several band widths to allow for wave length differentials In the 
absprptlon scattering characteristics of the slab and for band-dependent values of the Inci- 
dent flux. The slab may also be modeled with several different layers of materials to study 
the reflective behavior of a composite slab. This program uses Gaussian quadrature for the 
approximation of the scattering integrals and the overall reflectance of the slab. The program 
computes the beam Intensities at equally spaced stations, ©modes, throughcMjt, the slab. The 
output contains the Iterated beam intensifies for all stations within the slab as well as the net 
radiative flux and net radiative flux divergence for all nodes. The program Is able to handle 

*A discussion of the theory, program equations, convergence and results are given 
m Radiative Transfer In Highly Scattering Materials - Numerical Solution and Evaluation 
of Approximate Analytic Solutions” by Kenneth C. Weston, et al. 



non-unity Index of refraction cases and either isotropic or non-isotroplc scattering. 

The best way to describe this program’s capabilities and limitations is to present 
its necessary input and resultant output. Three specitrl cases have been run and are 
discussed below as representative of the various output formats. Immediately below is 
a description of all input data including their programmed range of values (the order of 
presentotion corresponds to the instructions found on the first page of the program listing): 
NDS - One of the most crucial decisions when setting up the data is the number of nodes 
to use for the slab. In an earlier paper (see footnote on page 1), the minimum number of 
nodes was found to be dependent on the optical thickness of the slab, its albedo, and th ‘ 
quadrature used; 

NDS>l+To(^^^^i^)forall!. (A1) 

Ui* 

As an example, with unity index of refraction (this program uses 16fh quadrature when 
the index is one), [u= .9995, and 5, a minimum of 49 nodes would be necessary to 
obtain a smooth convergent solution. It has been our experience, however, that accuracy 
of the solution increases as the number of nodes is increased. It is suggested to use 2 to 3 
times the minimum number from (Al) to insure high accuracy. Range: 11 751 In incre> '.i‘.nts 

of 10. 

THICK - The depth of the slab is arbitrary as long as the necessary optical thickness is 
maintained. Constant values of the scattering and absorption coefficients may be used for 
a variety of different optical depths with the only necessary adjustment to the data being in 
the thickness of the slab. For a multi-layer analysis, the slab is divided into five layers, 
each having a depth of 1/5 of THICK . Range: greater than zero . 

TOL - The program iterates back and forth through the slab solving for new values of the 
intensity field until the sum of the squores of the differences between the old and new 



values (PARAM ) is less than the set tolerance ( TOL ). For no changes in the fifth signi- 
ficant digit in the field a tolerance of 10 is necessary. Range; greater than zero, 

test - The program will abort if it reaches TEST iterations before obtaining a solution . 
Through some careless input of data a non-convergent situation may be encountered 
which would lead to an infinite loop. TEST terminates computations in this event. This 
variable may also be set to a very small number to determine initially if the input data 
Is in need of adjustment before making a long, expensive run. Range: greater than zero, 

PRINT - For the standard case the output contains information only from the last iferatic,; . 
One of the bands may be designated as a special band which informs the computer to print 
the information for the first, second, every PRINT, and the final iteration. For example, 
if PRINT = 10 and the computed PARAM became less than TOL on the 53rd iteration, the 
program would print out all data on the following iterations: 1, 2, 10, 20, 30, 40, 50, 
and 53. Range: Greater than zero . 

NONDM - Dimensionless data may be requested and will only appear for the final 
iteration , Beam intensities are divided by the Incident beam intensity . Radiative flux 
is divided by the incident flux. Flux divergence is divided by the absorption coefficient- 
incident flux product. In the case of different values for the absorption coefficient through 
the slab, the dimensionless values for the flux divergence will be based on the coefficient 
for the 5th, or bottom, layer. Range: Oor 1. 

BNDS - This variable fells the computer which bands to run and which one is a special 
band. If only a single standard band is being run, this card would contain a ”1” in the 
first column. If the first, second, and fourth bands were to be run with the second band 

as special, the data would be "1201". This data must be left justified on the card. Range: 
0, 1, or 2 for each of 10 bands. 



TEMPD - This sets the temperature distribution through the slab. Four options are avail- 
able: 1) The temperature at the front wall is used through the slab with the back wall 
temp independently set; 2) The back wall temp is used throughout with the front wall 
independently set; 3) The temperature falls or rises linearly from the front to the rear 
surface temp; and 4) The temperature for each node is individually read in. The first 
three options require reading in only the front and rear surface temperature while the 
fourth requires temperatures for all nodes. Range: Integers one through four. 

ISOT - The scattering function is set by this variable, (See paper listed in page 1 foot- 
note) Maximum forward scattering is achieved by a value of +1 while max backwards U 
-T. Isotropic scattering occurs when this variable is 0. “K " on the same card determines 
whether this variable is read in for each layer of the slab or whether the first value is 
repeated for all 5 layers. Range: -1 to 'H. 

M - The index of refraction outside the front boundary is assumed to be 1 . The slab Itself 
may take on a different index which may be band-dependent. Range: Greater or equal 
to 1 . 0 . I 

RB - This data may be used for either unifonm or Fresnel back surface reflection . The value 
of K determines which back surface condition is used and whether one value is used for the 
slab or ten band-dependent values are read in . If the Fresnel relationship is used, the 
real component (Nl)and the imaginary component ( Kl) of the substrate refractive index 
must both be read in . The program also checks to see if the Fresnel components are within 
the range prescribed for the approximation the program utilizes in this case, and will flag 
the output if the data is out of range. Range; For 0-> 1; for ^ and positive. 

QO “The band-dependent incident flux may be set for any value greater than zero. 



LAMBDA . - Tliis is the upper wave length limit in centimeters of a particular band. The 
program finds the lower limit to the band from the data on the previous band (for the first 
band it assumes a lower limit of zero). The width of the band determines the amount of 
energy introduced into the slab by emission in that particular band. Range: 

■ ABSCO ~ The absorption coefficient may be any non -negative value. 

SCATCO The scattering coefficient may also be any non-negative value. As described 
in the program listing, a set of the last few cards must be read in for each band up to the 
largest band num'ber referred to in BNDS . For example, if BNDS contains "101201'^ 
information must be read in for all of the first six bands. 

The program listing contains further explanation on the format of the input data 
cards and their order. Unless otherwise stated one card must be used for each card number 
referred to in the listing. 8A end B must be separate cards from 8. The data on the first 
six cards may be punched in any format but should be left justified. 

The final form of the program is presented below. This has been run on a XEROX 
SIGMA 5 with FORTRAN IV. Necessary adjustments must be made to the input/output 
statements before the program may be used on another system. Comment cards have been 
inserted strategically in the program to assist the alert and courageous user attempting to 
understand the logic behind the statements. 

Three test cases are also presented below as an aid to understanding the various 
output formats. The first case involves two bands with different indices of refraction. In 
this particular case the data for two cases was Identical except for the index of refraction 
and as such were combined Into one slab with the individual band output studied. Looking 
at the output, the data common to all bands is printed out first along with the number of 



the different bands to be run. The temperature for each node is printed next as a double 
check for the programmer. The scattering characteristic JSOJ is printed for each layer 
followed by the units for the beam intensifies, flux, and flux divergence . 

The next page of the output starts with the data peculiar to Band 1 . The absorption 
(K) and scattering (S) coefficients ore printed by layer. The Gaussian directions and 
weights for the eight forward directed beams are printed out since they are a function of 
the index of refraction (because the program incorporates discontinuity error avoidance). 
.With non-unity index of refraction three quadrature formulae are used to avoid discontinu- 
ity errors (as discussed in Appendix B of the Semi-Annual Status report for the period 
January to June 1974). 

The wave length interval is computed as described above. The input data is read 
in centimeters and internally converted to microns. 

F(0 - IT) describes the black body fraction of the emission at each node for this 
band. Since this fraction is dependent on temperature and band width, all the nodes are 
printed. In the case of zero temperature, garbage may appear in this output although this 
in no way affects the accuracy of the calculation of the intensify field. 

CRT is the cosine of the critical angle for total infernal reflection at the front 
boundary while CRTDG is CRT converted to degrees, 

RFll is the reflectance on the inside surface for each of the eight Gaussian direc- 
tions as printed out above, RFIO is the reflectance on the outside of the slab as computed 
for each direction corresponding to a single 16th order quadrature application. The reason 
for the difference between the quadrature formulae directions for RFLQ and RFLl will be 
discussed below. 

Information on the final iteration is printed next. Since this is a standard case, 
only every 10th node is printed. From 1 to 101 are listed the intensities of the forward 



directed beams and from 101 to 1, those for the backward directed beams. 1 1 corresponds 
to the beam intensity in the first Gaussian direction . PARAM is printed to double check 
that the solution did converge to the desired accuracy. 

The dimensionless data is also printed for only every 10th node. This calculation 
is followed by a reflectance field calculation/ with the following explanation* The 
Gaussian quadrature combinotion printed at the beginning of this band output was calcu~ 
loted so as to avoid a large error associated with the critical angle of total internal 
reflection. The reflectance field calculation (only used with non-unity index of refrac- 
tion) uses the 16 internal directions corresponding to the 16 external directions in 16th 
order quadrature. The program iterates across these new directions within the already 
established intensity distribution. This is done to calculate, by using RPLO , the value of 
the reflected beams on the outside of the slab in those directions corresponding to 16th 
order quadroture so that the intensities may be readily integroted into a very accurate 
estimate of the overall slab reflectance. Immediately following this print out is the 
dimensionless value of the intensity field on the outside of the slab in the eight backward 
directions. From this the band -dependent slab reflectance is calculated and printed. This 
reflectance value is saved from each band to form an overall reflectance for the slab. 

On the next pages of the output the data for the second band, which follows the 

formot described above for the first band,is given. Note that the Gaussian directions and 

♦ 

weights are not the same for the two bands. 

The second special case was run to test the composite data from a multi-band slab 

analysis. The Incident flux for each band was chosen such that the composite data corre- 

2 

spends to a single black band case with 459.646 watts/cm"^ incident flux. A comparison 
was made with a previously run single black band case and the two were found to match 


well wUhin allowable errors. Also in this test case the second band was run as a special 
band to demonstrate the output format for a special band. The output starts with the 
speciol band information. Note here that the black body fractions are non-zero due to 
the non-zero temperature distribution. The printout for the special case is similar to the 
standord except it contains the calculated reflectances across the front boundary from the 
outside in (RO) and the inside out (Rl). Also printed is the starting routine where the 
Kubelka-Munk approximation is used to provide a first value for the intensity at each 
node from which to iterate the intensity field. The printout contains the information on 
all nodes for the first, second, every 7th, and final iteration, although most of these 
pages were left out here for brevity in the presentation . 

On the following pages are the output for the first, second and third bands. Notice 
that the condensed information from the second, or special, band is repeated in its appro- 
priate place in the output. After the third band the overall reflectance Is printed out 
followed by the summed net radiative flux and net radiative flux divergence from ail three 
bands. 

In the third test case a slob is divided into two layers os described on the printout. 

This single band case is njn as a special band so that the intensity distribution may be crossly 

studied node by node near the interface between the layers. The dimensionless value for the 

final iteration is presented in its entirety since this particular case generates a most inter- 
/ * 

esting intensity field. 

.As a role much more detailed information is printed out on a special band than on a 
standard band but the computer time is proportionately longer since the printer is working 
longer. 

The program has been diligently checked out and has shown no visible errors within 
the limits of the cases studied. 
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COMMON Fj IPOS/ INEG# IP eSPi INEQP/SAVE . y BDOOOl 

real A|8I#AA, AliA2/BB/BliB2#BET^Cl#C2iC3#CRT>CRTDG/D0R(75i l/OQROY^ BDOOt 
IDQRSAV! 751 )iDYiDIV^ DMLQR j DMLDQR/ F ( 10# 751 » # F I# I P6S ( 8 # 75 1 1/ IPOSN# BD002 

21NEG( 8#751 )# IPOSP ( 3# 751 I i INEGP ( 8# 751 » # KI UO)#KF( 10# 5 ) # KKM/ LAMBDA ( BD003 

31O)#1S0TI5)#MU(8) jMUP(8)#MUe#MUU#Nao)#OARFL#0ARFLl#0ARFU2#PARAM# BD004 
4PHASE(8#8)#PHASC0iPHASKa/Q(75l )#QT#ORR#QTP#ORP#OOUO 1#Q0O#ORSAV(75 BD005 
51 l#RB( 10)#RBB/RIN/RFLI ( 8) #RFL0( 8 »#RJ#RI#R( iO)#R0#S#Sl#SF U0#5 l#Sia S0n06 
real save (308# 11 )#SUM#SlNMU#SINMUe#SKM#T(751)#TBW#TFW#TK#TeRM# BD007 

1THICK#T0L#V# VALUE# Y # A3# A4 # C4# FRFLBUO I # FRFL I (10# 8 > # FRFLO ( 8 » # 00008 

2DML(8)#KFF#SFF#NN#RF(8)#MU1(8)#AU1(8| BDOOSl 


INTEGER Ae0RT#BND(ll)#BNDS(10)#I#IJ#II#IK#It#lN,llI#IJI#J#JL#JK# BD009 
IJM# JP# JS# JT# J2# JR# JQ# JJJ#K#KK#KJ#fiaga*^f'R£8'^'l-L#LJ#LLL#LnT(6l#Mt1# BOO 10 

pmm nds#nondm#print#prtvar#puacE(8)#tempd#test#ja soon 


DATA INPUT# BY CARD 

1. NOS - NUMBER OF N0DES 

2. THICK - THICKNESS IN CM 

3. TOL - PARAM TOLERANCE 

4* TEST - MAXIMUM NUMBER OF ITERATIONS 

5, PRINT - FOR SPECIAL BAND/ PRINTS FIRST# SECOND# EVERY (PRINT)# 

AND LAST ITERATION 

6« N6N0M - O-N0# 1-YES 

7* 0NOS - TO SET UP BANDS# lOIli 0-DBN»T WORK# 1-STANDARD# 2-SPECIAL 


8* TEMPO - temperature DISTRIBUTION# i-FW THROUGHOUT# BW| 

2-FW# BW THROUGHOUT! 3*ipW LINEAR 0Hj 4-READ IN DlSTRlBUTIONt 
IF 1*3, ,1 

8A* TEMPERATURE FRONT WALL 
8B< TEMPERATURE BACK WALL 
IF 4, 

8A - 8? FRONT WALL TO BACK# BY NODE# lOF8t3 


9, K# IS0T ( 1 1# 5F10* 7 ) - ANISOTROPY BY DEPTH. K-0# FIRST VALUE REPEATED# 

K-l# ALL 5 VALUES ON CARD* 

10* K#N( M#10F7.3) - INDEX OF REFRACTION OF MEDIUM BY BAND* K AS ABOVE* 

11* K#RSm# 10F7.3) - RB# OR NI I F USING FRESNEL RELAT I ONSH IPS# BOTH BY BANO* 
» « T'-.n. n T*- MR R- FR'^^'Nn V‘l HP®* 



-- * 



^ <Ji '>1’ <v>.r LAIkU. 1 f'l r AUL. iU vHUUi^w 


50'“. 

C 


K - 2/ NI AND KI REPEATED# ALTHOUGH KI 0N SECOND CARD (10F7«3) 

• 

51* 

c 


K - 3# NI READ IN# KI REPEATED. K • 4# NI REPEATED# Kl READ IN* 

52. 

c 


K •- 5# NI and KI READ IN. . 


53. 

c 




54. 

c 

note; cards 12/ 13^ AND 13A (IF REQUIRED) MUST BE REPEATED FOR EACH BAND(J). 

55* 

c 




56. 

c 

,12. 

K#00( J)#LAMBDA( J) ( I1#F15.7#E10.5) - INCIDENT RADIATIVE FLUX AND UPPER 

57. 

c 


LIMIT OF BAND. 


58. 

c 




59. 

c 

13. 

KF(JM5F10.6) • ABSCe#SCATC0. USING THE K FROM CARD 12# 


60* 

c 


K-0# 0NE VALUE F0R EACH ON THIS CARDi K-l# ONE VALUE FOR ABSC0 

# 

61 . 

c 


5 VALUES FOR SCATCO 0N NEXT CARO ( 5F10. 6 ) | K-2# 5 VALUES FOR ABSC0#, 

62. 

c 


one VALUE FOR SCATCO ON SECOND CAROl K»3# 5 VALUES FOR EACH. 


63 . 

c 




64. 

c 


' 


65. 

c 

14. 

LAST CARD - MUST HAVE '4» IN FIRST COLUMN OF FIRST CARD# FOLLOWED 

BY 

66 « 

c 


BLANK CARD. 


67. 

c 




68. 

c 




69. 



DATA MUU ) #MU( 2) #MU(3) #MU( 4J/.09501250984# .28160355078# .4580167776 

BD012 

70. 



16# .61787624440/ 

6D013 

71. 



DATA MU (5 )#MU( 6) #MU( 7 )#MU( 8)/. 75540440836#. 865631 20239# .9445750230 

BD014 

72. 



17# .98940093499/ 

BD015 

73. 



DATA A( 1 ) # A ( 2)# A( 3) # A( 4 )/.l 8945061046# *18260341504# <16915651940# 

BD016 

74. 



1 .14959598882/ 

BD017 

75. 



DATA A(5)#A(6)#A(7)#A(8)/.1 24628971 26#. 0951585 1168# *06225352394# 

BD018 

76. 



1.02715245941/ 

BD019 

77. 



S » 5.6699E-12 

0DO2O 

78. 



D0 4151 = 1#8 

600201 

79. 



MUl { I) » MU( I) 

B00202 

80. 


415 

AUl ( I) » AU » 

BD0203 

81. 



0ARFL1 « 0* . 

BD030 

82. 



0ARFL2 * 0. 

BD031 

83. 



ABORT « 0 

BD032 

84. 



INPUT NOS# THICK# T0L# TEST# PRINT# NONDM 

BD033 

85. 



OUTPUT NDS#T0L# TEST# PRINT# N6N0M 

60034 

86* 



PRINT 901#THICK 

BD0341 

87. 


901 

FORMATl IX# 'THICK »'#F8*5#' CM*) ... . 

BD0342 

38. 



DO 4001 ■ l.#NDS 

BD035 

89. 



QRSAV(I) * 0. 

. BD036 

90. 


400 

DQRSAV( I ) « 0. 

BD037 

91 . 

c 




92. 

c 

SET UP BANDS. 


93. 

c 




94. 



READ IjBNOS 

BD038 

95. 


1 

F0RMAT(1OI1) 

60039 

96. 



IJ « 1 

BD040 

97. 



BNO 1 1 ) ■ 0 

B0041 

98. 



DO 21 * 2# 11 

BD042 




’ ' r f T 1 1 - n • 

pnn- ?i 


1 uu • 



wi'lU \ i 1 U 


U U J T 

ion 



IF(BNDS(I » 1) - 1)2/3j4 


BD044 

102« 


4 

BND( 1 ) - I - 1 


BD045 

103* 


3 

BNOl I J + 1 ) » I - 1 


B0046 

104* 



IJ = IJ + 1 


BD047 

105* 


2 

CONTINUE 


BD048 

106* 



II • 2 


BD049 

107* 



IF(BND(1> )7/7i6 


BD050 

108* 


6 

11 = 1 


BD051 

109* 



IF! IJ - 2)5^5i7 


B0052 

no- 


5 

IJ « 1 


BD053 

un 



BND(2) « 0 


B0054 

lia- 


7 

PRINT 8/BNDSnBN0( I )/I » lilJ) 


BD055 

na* 


8 

FORMATt' BNDS '/10H/» SPECIAL — 

ni2/» STANDARD --*/10I3/l 

B0056 

114- 

C 


• 



115- 

c 

SET UP temperature DISTRIBUTION 

t 


116* 

c 





117- 



INPUT TEMPO 


BD057 

118- 



OUTPUT TEMPO 


BD058 

119- 



iFfTEMPD " 3l9i9>lQ 


BD059 

120- 


9 

INPUT TFWiTBW 


BD060 

121- 



T( 1 > “ TFW 


BD061 

122- 



T(NOS) .« TBW 


B0062 

123- 


- 

IF(TEMPD 2)11/12/13 


BD063 

124- 


11 

TK « TFW 


B0064 

125- 



GO TO 14 


SD065 

126- 


12 

TK w TBW - .. 


BDQ66 

127- 


14 

DO 151 » a^NOS-1 


BD067 

128* 


15 

T( I ) • TK 


BD068 

129- 



GO TO 16 


BD069 

130- 


13 

TK e {TFW - TBW)/. (NOS - 1- 1 


B0070 

131- 



DO 171 « 2/NDS - 1 


BD071 

132- 


17 

T( I ) « T( I - 1 ) - TK 


BD072 

133- 



GO TO 16 


BD073 

134* 


10 

READ ISnTI I )/I “ 1/NDS) 


BD074 

135* 


18 

F0RMAT( 1 10F8.3) ) 


BD075 

136- 


16 

PRINT 19/ ( T ( 1 )/ I « 1/NDS ) 


8D076 

137* 


19 

FORMATi* TEMP » '/ 10F9 - 3/ 1 X/ M KELVIN » 

»/t8X/10F9.3/l ) 

BD077 

138- 

c 



j 


139* 

c 

SET 

UP ANISOTROPY DISTRIBUTION 

1 


140* 

c 





141- 



READ 20/K/ IS0T 


BD078 

142- 


20 

FORMATt I1/5F10.7) 


BD079 

143- 



IF(K)2l/2i/22 


BD080 

144* 


21 

DO 231 » 2/5 


BD081 

145- 


23 

ISOTd) « IS0TI1) 


0DO82 

146. 


22 

PRINT 24/ISeT 


BD083 

147- 


24 

FORMATt ’ IS0T « »/5F11.7/) 


B0084 

148- 

c 





149. 

c 

SET 

UP INDEX OF REFRACTION 





l^X * 



U w 

152* 


26 F0RMAT( Ilil0F7*3) 

BD0S6 

153« 


IF(K)27>?7/28 

BD087 



27 D0 291 « Z,i0 

BD08S 

155* 


29 N( n N( 1 ) 

BD089 

156* 

C 



157* 

c 

SET UP RB 


158* 

c 



159* 


28 read 32/K/RB 

B0090 

160* 


32 FORMAT ( I1>10F7*3) 

BD091 

161 • 


KFRES « K 

BD0919 

162* 


IF(K - 1}33j34^36 

BD0920 

163* 


36 READ 37/KI 

BD0921 

164* 


37 FaRMATl 10F7*3) 

B00922 

165. 


IF(K - 5)30#42/42 

BD0922 

166* 


30 IF(K - 3)3Bi39/38 i 

BD0923 

167* 


38 D8 401 s 2il0 

BD0924 

163* 


40 R8S I) = R3{ 1 J 

BD0925 

169* 


IF(K • 4>39/42i39 

BD0926 

170* 


39 D0 411 « 2#10 

000927 

171 * 


41 KI ( I 1 « KI ( 1 ) 

BD0928 

172* 


42 PRINT 43 

BD0929 

173* 


43 F0RMAT11X//' FRESNEU"S CRITERIA*//* BAND VALUE*/» 

BD0930 

174* 


00 441 « 1/10 

BD0931 

175* 


IFINI I ) )44/ 44|63 

BD0932 

176* 


63 IFIRBd) )44/ 44i64 

0DO933 

177* 


64 VALUE « (Re(I)¥RB(I) + K I ( II ¥KI ( I M *5/N ( I ) 

BD0934 

178* 


IFIVALUE - 3*3)72/72/73 

BD0935 

179* 


72 PRINT 74/ 1/ value 

B00936 

180* 


74 F8RMAT(2X/ I2/3X/F7.2/ ' NOTES OUT OF RANqE OF FRESNEL APPROXiriATI 

BD0937 

181* 


iON’ ) 

BD0937 

182* 


GO TO 44 

800938 

183* 


73 PRINT 79/ I/VALUE 

BD0939 

184* 


79 F0RMATI2X/I2/3X/F7.2) 

BD0940 

185* 


44 CONTINUE 

BD0941 

186* 


GO TO 34 

BD0952 

187* 


33 00 351 » 2/10 

BD0953 

188* 


35 RQ( I ) • R0( 1 ) 

BD0954 

189* 

c 

1 


190* 

c 

SET UP BAND DATA 


191 * 
192* 

c 

34 J = 1 

BD0955 

193* 


62 READ 45,KiO0( J)/LAMBDAI J)/ tKFIJ/I)/I ■ 1/5) 

B0096 

194* 


45 FORMAT ( 11 / F 1 5 * 7/ E 1 0 • 5/5F 1 0 • 6 1 

BD097 

195* 


IF(K •* 4)60/61/61 

BD098 

196* 


60 IF(K •* 1 )46/47#47 

BD099 

197* 


46 SF( J/1 ) = KF( J/2) 

BDIOO 

198* 


DO 491 • 2j5 

BDlOl 

199. 


SF( J/ I } « SF( J/ 1) 

80102 

200* 


49 KF( J/ I ) * KF (J/ 1 ) 

BD103 

'*■ T 1 • 


rn TPi "n • 

ani n 


i # u 1 


I%h. U W/iWi 


w U I U 


203". 


48 

F0RMAT(5F1O.6) 

BD106 

E0^< ' 



IF{K - 2)51/52/50 

BD107 

205* 


51 

DO 531 » 2/5 

00108 

206. 


53 

KF( J/ I 1 « KF( J/1) 

. 80109 

207* 



G0 TO 50 

BDllO 

208* 


52 

DO 541 « 2/5 

BDin 

209. 


54 

SF(J^ I ) « SF( Jil) 

BD112 

210. 


50 

J-J + 1 

B0113 

211 . 



IF(J - 11162/61^61 

00114 

212* 

C 




213* 

C 

SET UP BOUNDARIES FOR NODE SECTIONS 


214. 

C 




215. 


61 

KK • NDS/5 

BD115 

216* 



LMT12) » KK , 

60116 

217. 



LMT( 1 ) - 0 i 

80117 

218. 



LMT(3) = LMTJ2) + KK , 

B0U8 

219. 



LMT(4) « LMT(3) + KK + 1 

BD119 

220. ' 



LMT(5) * LMT(4) + KK 

BD120 

221 . 



LMT(6) = LMT(5) + KK 

BD121 

222. 

C 




223. 

C 

SET UP BLACK BODY FRACTIONS ACCORDING TO BAND WIDTH AND TEMPERATURE 


224. 

C 



1 

225. 



IJI « BNO( IJJ 

BD122 

226. 



III » BN0^2) - 1 

BD123 

227. 



IF( III • 1 )610/611/611 

BD124 

228. 


610 

III » 1 . 

B0125 

229. 


611 

DO 102MM a IJI/III/-1 

B0126 

230. 



FI a LAMBDA(MM) 

BD127 

231. 



00 102LL - 1/NDS 

B0131 

232. 



IFIT(LL) 1102/102/25 

B0131 1 

233. 


25 

V « 1.43879/F1/T(LL) 

BD132 

234. 



IFI V - 2. )103/104j104 

B0133 

235. 


103 

F(MM/LL| a 1. - 15./(3.141593|»^44»V¥»3^(1./3. - V/8. + V^^V/60# - V 

BD134 

236. 



l¥¥4/5040. + V^f*6/272160. - V-^» + 8/13305600. ) 

B0135 

237. 



G3 TO 201 

BD136 

238. 


104 

SUM a 0. 

BD137 

239. 



DO 105LLL a 1/5 

B0138 

240. 


105 

SUM a SUM + EXP(-LLL^V)/LUL'f^4^n (LLt^V + 3*»¥LLL^V ♦ 6*)¥|^UL*V + 

BD139 

241. 



16. ) 

BD140 

242. 



F(MM/LL1 a SUM»15./{3.141593)^‘¥4 

BD141 

243. 


201 

IF(MM - IjniOO/102/102 

BD142 

244. 


100 

F(MM + 1/LL) » F(MM + 1/LD - F(MM/LL) 

BD143 

245. 


102 

CONTINUE 

BD144 

246. 



DY a THICK/INDS - 1- ) 

BD145 

247. 



PRINT 902 

BD1451 

248. 


902 

F0RMAT(lX//» INTENSITIES IN WATTS/CM»»2/STERA0I AN ' // ' FLUXES IN WA 

BD1452 

249. 



1TTS/CM+¥2'//» FLUX DIVERGENCE IN WATTS/CM^*3 » 1 

BD1453 

250. 

C 




251. 

c 

MAIN DO-LflOP FOR RUNNING THROUGH EACH BAND 





r 


U'.J I '■ 


1 1 i 1 o 


w 1 


u -J * 

254% 

255 * 

256 * 

257 * 

258 * 

259 . 

260 

261 

262 

263 

264 

265 

266 
267 
263 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 ' 
302 ' 


C 

c 

c 


c 

c 

c 


SET UP variables PECULIAR TO EACH BAND 


KK * BNOI IK ) 

NN = N{KK|¥N(KK) 

DO 6801 « 1#B 
MU ( I 1 « HUl ( I ) 

680 A ( I ) » AUl ( I ) 

R8B » RB(KK) 

IF(NN -- 1 )850/850#860 
860 RMUC » (1* - l./NN)^^t,5 
TERMl = ( 1* - RMUC)».5 
TERM2 « (1* + RMUC)¥*5 

CALL BRKDWN{RMUC>TERM1j»TERM2#AjMU) _ 

850 IF(KFRES - I)113>ll3jn4 
114 I «> KK 

FRFLB( I ) « 0* 

IF(N( I ) )83i83i84 
84 IF( RB( 11 )83f 83^101 
101 DO 219J « liS 

A1 » 1 RBI I )^MU( J ) - N( I) )¥¥2 
A2 » (RBUI^MUIJ) + N(I))4¥2 
A3 - 1 R6 1 1 ) • N( I I ¥MU( J) )¥¥£ , 

A4 = (R0( I) + N( I )¥MU( J1 1¥¥2 
Bl » K1(I)¥KI(I) 

92 » B1¥MU( J1¥MU( J> 

FRFlKI/J) “ *5¥((A1 + B2)/(A2 ♦ B2) ♦ U3 + B11/U4 + BD) 
219 FRFLBil) » FRFLB(I) + A ( J i ¥MU i J ) ¥FRFLI ( I # J ) 

FRFLBU) * FRFLBt I )/*50l5l552 
83 RBB « FRFLB(KK) 

113 Q00 « Q0(KK) 

IPeSN » O0O/3»141593 
IFUK - l)66i66>67 


PRINT HEADING AND INFORMATION 


66 PRINT 68#KK 

68 F0RMAT('l SPECIAL BAND « BAN0*#I2//) I 

Q0 TO 69 

67 PRINT 70^KK 

70 FORMAT! »1 BAND'rlE//! 

500 FORMAT! lXy//» ITERATION »/I3//) 

69 PRINT 7liN(KK)iRBBiQ00i tKF(KK/JI# J - 1 # 5 1 i ! SF i KK^ J ) / J - 1^5) 

PRINT 853/ (MU(J)jA(JI^ J • 1/81 

853 F0RMAT!' GAUSSIAN I NTEQR AL • /4X/ 4HMU » S/ 7X# » WE I GHTS » /8 ( 1 X/ FIO • 8/ 2X/ 


71 


1F10*8/M 

FORMAT! » INDEX -'/F8.3/’ RB «»/F6«3/» QO »»/F16.7/» WATTS/CM*¥2 « / 
1' K «»j5F11*6/» CM»^-1V* S "»/5F11«6/* CM»¥-1'//1 
IF!KK - 1 1904/904/905 


BD1461 

BD1462 

BD1463 

BD1464 

BD1465 

BD1470 

BD1471 

BD1472 

BD1473 

BD1474 

BD1475 

BD14S4 

BD1491 

BD1492 

BD1493 

BD1494 

BD1495 

B01496 

0D1497 

BD1498 

BD1499 

BD1500 

BD1501 

BD1502 

BD1503 

BD1504 

BD1505 

BD1506 

BD151 

BD152 


BD153 

BD154 

B01.55 

BD156 

BD157 

BD15S 

BD159 

BD1591 

BD1592 

B01593 

B0160 

BD161 

BD1611 

°ni ' \ o 






Uu lu 9 u6 

>-'U ], l. J 

305 

904 

A1 ** 0* 

B01614 

306* 

906 

A2 * LAMBOA(KK) 

B01615 

307- 


A3 » Ai^lOOOO* 

BD1616 

308* 


A4 a A2-V10000* 

BD1617 

309- 


PRINT 903^ Al/ A2i A3i A4 

BD1618 

310« 

903 

FOR^ATtix//’ Wavelength interval; *iEio*5#* - »jE 10>5/» cn'/23Xi 

BD1619 

311* 


IEIO.5^' - '/El0.5#t HICRONSV) 

6D162 

312. 


PRINT 701MF(KKi J)/J • liNDS> 

6D1621 

313. 

701 

FORMATdX,' F(0 - LT ) ; V( 1X> 10EXl*5/n 

BD1622 

314* 


IF( IK - 1 )75/75/76 

B01623 

315* 

76 

IF(KK - BNDI 1) )75>78#75 

BD163 

316. 

317* 

C 

C IF 

ALREADY RUN AS SPECIAL BAND> PRINT SAVED NUMBERS 


318* 

C 



319. 

78 

PRINT 82 1 

B0165 

320* 


PRINT 80 

BD1651 

321* 

80 

FORMAT ( ' NODE d6X# ' 1 1 * i 9X/ ' 1 2 * / 9Xi * 13 * / 9Xi * 1 4 d 9X# ' 15 d 9X# * 1 6 ' i 9Xi 

BD166 

322. 


l'l7d9Xi *I8»i9Xi 'QRW8X/ ’DQRDY*/) 

00167 

323. 


PRINT 81/ ( (SAVE{IiJ)^J - 1/11)/I • XfL) 

B0168 

324* 

81 

FORMAT! ( 1X/F4.0/2X/9|E10«5/ 1X)/E11.51) 

60169 

325. 


PRINT 82 

BD170 

326. 

82 

FORMATdX//) 

8D171 

327- 


PRINT 81/ ( (SAVEi I/JJ/ J « I/IU/I « L+l^L + U) 

BD172 

328. 


OUTPUT R(KK) 

801721 

329* 


IF!N0NDM»65/65/S5 

BD173 

330. 

85 

PRINT 86 

BD174 

331. 

86 

FORMAT! IX///’ DIMENSIONLESS'//) 

BD175 

332. 


PRINT 80 

B0176 

333. 


PRINT 81/ ( (SAVEil/ J)/J 1/llJ/I - L + L+l/L + L + L) 

BD177 

334. 


PRINT 82 

B0178 

335. 


PRINT 81/ ( (SAVElI/ J)/ J « l/ll)/I » L+t+L + l/U + L + L+L) 

BD179 

336. 


GO TO 65 

BD180 

337. 

C 



338* 

C calculate RFLECTANCE for inside GAUSSIAN ANGLES 


339. 

C 



340. 

75 

RIN = 1*/N(KK) 

BD181 

341 . 


CRT =* ARSIN(RIN) * 

BD182 

342. 


CRTDG = CRT»l80*/3. 141593 j 

BD183 

343. 


OUTPUT CRT/CRTDG 

B01831 

344* 


DIV • COS(CRT) 

B0184 

345 1 


DO 1101 =1/8 

BD185 

346. 


IF(MU!I ) - DIV)ill/lll/112 

BD186 

347. 

112 

SINMU » (1* - MU( I )¥MU( I 

BD187 

348* 


SINMU0 « SINMU¥N(KK) 

BD188 

349. 


MUO « !1* • SINMU0¥SINMUO)¥¥.5 

BD189 

350. 


A1 = SlNMUO¥MU(n + MU0¥SINMU 

BDian 

351 • 


A2 « SINMU0¥MU(I» - MU0#SINMU 

BD191 

352. 


B1 • MU0¥MU(IJ - SINMU¥SINMU0 

B0192 

353. 


B2 « MU0¥MU(I) + SINMU^SINMUO 

BD193 

ri ^5 ' • 



RDI QA 


356 


111 

RFLH I ) = 1. 


B0196 

357. 


no 

CONTINUE 


B0197 

358. 



PRINT 55#RPLI 


BD1971 

359. 


55 

FORNAT(' RFLI « 'i8F10.7/» 


BD1972 

360. 

C 





361. 

C 

CALCULATE OUTSIDE QAUSSIAN ANqLES AND C8RRESP0N0INQ REFLECTANCES 


363. 

c 





363. 



D0 1151 » liS 


B019S 

364. 



SINNUe « 11. - MUl ( I )4MU1 ( I ) M+.5 


BD199 

365. 



SINMU .. SINMUB/N(KK) 


BD200 

366 . 



MUP(I) » (1* • SINMU¥,SINNU)¥^.5 . 


B0201 

367. 



A2 * SINMU0#NUP(I) - MUnn^SINMU 


BD202 

368. 



A1 » SINMU04MUPI I) + MUnn^SINMU 


B0203 

369. 



B1 » MUim^MUP(I) - SINNU^'SINMUe 


BD204 

370. 



82 « MUnn^^MUPlI) 4 SINMU^SlNMUe 

1 

B0205 

371. 


115 

RFLS(I) « .S^tAE^Aa/Al/AD^^d. + Bl*8l/B2/B2 J 


BD206 

373. 



PRINT 56iRFL0 


B02060 

373. 


56 

F8RMATI ' RFUe « d8F10*7/I 


BD2060 

374. 

c 





375. 

c 

KUBELKA MUNK STARTINQ ROUTINE 



376. 

c 





377. 



IF(KFRES - 1)106/106/108 


B02062 

378. 

\ 

loa 

DO 121JN « 1/8 ... 


BD2063 

379. 



A1 » (RB(KK)#MUP( JN) .• N(KK))»¥2 


BD2064 

380* 



A2 « (R0(KK).^NUP( JMI + N(KKM*Jf2 


BD2065 

381. 



A3 » IRBtKK) - N(KKl¥MUP(JMn*»2 


BD2065 

382. 



A4 * (RB(KK) + N(KK)4MUP( JM) )»42 


BD2065 

383. 



81 » KI(KK)»KI(KK) 


BD2066 

384. 



B2 a Bl¥NUP{ JM)¥MUPUM) . , . 


B02067 

385. 


121 

FRFlOUM) « .5iM(Al + B2)/(A2 4 B2» + (A3 ♦ Bl)/tA4 

+ Bi n 

BD2068 

386. 


106 

SI » 0. 


B0207 

387. 



DO 1201 a 1/8 . ... 


B0208 

3S8. 


120 

SI a SI + A( I )4MU( I )♦RFLI ( n 


BD209 

389. 



RI a SI/. 50.151552 


BD210 

390. 



R0 a 1. + NN4(RI » i. ) 


. BD211 

391. 


133 

Y a -DY 


B0212 

392. 



IF( IK - 1 ) 130/ 130/ 131 


B0213 

393. 


131 

IFUB0RT1125/ 125/ 13O 

1 

B0214 

394. 

c 



s 


395. 

c 

PR 

INT OUT starting ROUTINE IF SPECIAL CASE BR ON ABORT 

STATUS 


396. 

c 





397. 


130 

OUTPUT R0/RI 


BD215 

398. 



PRINT 132 . 


B0216 

399. 


132 

F0RMAT( 1X//3X/ tKUBELKA MUNK STARTINQ R8UTINE ' //2X, ' N0OE ' / 5X* ' OT ' , 

0D217 

400. 



110// ’0R»/9X/ 'OQRDY»/) 


BD218 

401 . 


125 

00 1221 a 1/5 


BD219 

402. 



III a LMT(I) + 1 


BD220 

403. 



JJJ a LMT( I + 1 ) 


BD221 

404 . 



SKM a SF(KK/I )'^.75 


BD222 

• nc; , 



'' 'M K'F f • T 1 ? 


Rn?5>3 


U U V. 


408* 

409* 

410. 

4U. 

4ia* 

413. 

414. 

415. 

416. 

417. 

418. 

419. 

420. 

421. 

422. 

423. 

424. 

425. 

426. 

427. 

428. 

429. 

430. 

431. 
432* 

433. 

434. 

435. 

436. 
437* 
433. 

439. 

440. 

441. 

442. 

443. 
444 * 

445. 

446. 

447. 

448. 

449. 

450. 
451* 

452. 

453. 
454 . 
455. 


^ 1 O i ! V ( t. . oi>n » I ^ 

BET = SIG/(KKM + g.^SKM) ! 

A1 » QO0*£XP(-SIQ^THICK)4U. - Re)*(BET*(l. + RBB) - 1. + RBB) 
A2 » 2.4 ( ( BET¥BET»J 1 . + RU^U. + RBB) + U* - Rl)*(.l* - RBB))* 
1S1NH(SIQ4THICK) + 2.*BET*(1. - RI *RBB ) *ceSH ( SI G^TH I CK I ) 

Bt « QQa*EXP(SIG*THICK)*( 1. - Re)*(BET*U. * RBB) + 1. - RBB) 

AA « A1/A2 
BB - BI/A2 
Oe 122J = IlliJJJ 
Y » Y + DY 

QT » AA4(1. - BET)4£XP(S1G4Y) + + BET ) 4EXP ( -SI Q*Y ) 

QRR = AA4(1. + 8ET)*EXP(SIG*Y) + - BET ) *EXP ( -S I Q4Y ) 

DQROY => -KKM4((3T + QRR I 
QTP » OT/3. 141593 

QRP « QRR/3. 141593 _ , 

DO 123JL « 1j8 1 

IPeS(JL/JI - QTP 

123 INEG(JLiJ) - QRP 

IF (IK - 1 ) 124/ 124/126 

126 IFtAB0RT)l22/122il24 

124 RJ a J 

PRINT 127/RJ/0T/QRR/0QRDY 

127 FORMATUX/FS.O/ 1X>2(E10.5/2X1/EU*5) 

122 continue 

' D3 2QOJK « 1/TEST 
PARAM a 0 . 

c . 

C FRONT SURFACE CALCULATIONS 
C 

DO 1501 - 1/8 
IPeSd/l) - INEGd/1) 

IF(MU(I ) - OIV) 150/150/151 

151 IPOS(I^l) » RFH ( I)*INEGI 1/ n + U* - RFLH I) l*IP6SN*NN 
150 CONTINUE 
C 

C MARCH TO BACK WALL 

C . . 

DO 1621 a 1/5 

00 161IL = 1/8 . 

DO 161IN =1/8 

161 PHASEIIL/IN) = 1. + ISOT { I )*MtM IL)*MU(IN) 

III » LMT( I ) + 1 
JJJ a LMT( I + 1) ■ 

IFlin - l)163U63d64 

163 III a 2 

164 SFF a SF ( KK/ I ) 

KFF a KF(KK/ I) ^ 
bet a SFF + KFF 
C2 a DY4SFF/2. 

C4 a DY4S4KFF/3. 1415934NN 

Off 1 ' I T T T . 1 1 I 


B0225 

BD226 

60227 

BD228 

BD229 

BD230 

8D231 

BD232 

BD233 

B0234 

BD235 

BD236 

BD237 

BD238 

BD239 

BD240 

BD241 

BD242 

BD243 

BD244 

BD2441 

BD245 

BD246 

BD247 

BD24S 


BD249 

B0250 

BD251 

B0252 

80253 


BD254 

BD255 

BD256 

BD257 

BD258 

B0259 

BD260 

60261 

BD262 

B0263 

BD264 

B02641 

BD2642 

an?" 6 



458J 



DO 162JL » li8 



BDP66 

459* 



tiUU, =* MU( JL) 



BD267 

460. 



SUM « 0. 



BD268 

461. 



Cl « 1. - BET'fOY/MuU 



BD269 

468* 



DO 165LJ « 1/8 



BD272 

463. 



PHASCO ■ PHASE(LJ/JL) 



BD273 

464* 



PHASK0 « 2. * PHASCO 



BD274 

465. 


165 

SUM « SUM + A(LJ)^(PHASCO¥lP0StLJ/J -1* + 

PHASKO^INESaJ/ J - 

n ) 

B0275 

466. 



TERM » IP0SUL/J) 



BD276 

467. 



IPOS<JL/J) » Cl¥-IP0S( JL/ J - 1). + (C2»8UM + 

C3)/MUU 


BD277 

468. 


162 

PARAM « PARAM + (TERM/IPOSN - IPOS « JL/ J 1 / IPOSN ) ♦♦a 


BD278 

469. 

C 






470. 

c 

BACK BOUNDARY CONDITIONS 




471. 

c 






472. 



IF(KFRES - 1 ) 128/128/129 


i 

BD2781 

473. 


129 

Cl » S + NNj^T(N0S)^^4/3«l4l593¥r(KK/J» 



BD2782 

474. 



DO 1161 « 1/8 



BD2783 

475. 


116 

INEG(IiNDS) = FRFLI IKK/ I UIP0S( I/NOSI + (1. 

r RBB)*C1 


B027R4 

476. 



GO TO 117 



BD2785 

477. 


128 

Cl « (1* - RBB)^S¥NN*T(NDS)4i»i4/3.1415934nF(KK/J) 


BD279 

478. 



DO 1701 « 1/8 



BD280 

479. 


170 

INEQ(l^NOS) « RBB¥lP0S( I/NDSl + Cl 



8D281 

480. 

c 






4S1. 

c 

NARCH TO FRONT WALL 




482. 

c 






483* 


117 

DO 1801 * 5/ 1/-1 



BD282 

484. 



DO 181IL « 1/8 . 



BD283 

485. 



DO 181IN » 1/8 



BD284 

486. 


181 

PHASEtIL/IN) » 1. - ISOT(II-»‘MUaU^‘MU(IN) 



BD285 

487* 



III = LMT( I + 1 I 



BD286 

488. 



JJJ » LMT( I) + 1 



BD287 

489. 



IFIIII NDS)182/ 183/ 183 



BD288 

490. 


183 

III » NDS - 1 



BD289 

491. 


182 

SFF « SFtKK/I ) 



BD290 

492. 



KFF » KF ( KK/ I ) 



BD291 

493. 



BET « SFF + KFF 



B0292 

494. 



C2 = DY»SFF/2. 



BD2921 

495. 



C4 = DY^S*KFF/3. 1415934NN 

1 


802922 

496. 



DO 180J * II1/JJJ/*1 

f 


BD293 

497. 



C3 » C4mT(J + 1)^»4*FIKK/J) 



BD2931 

498. 



DO ISOJL » 1/S 



BD294 

499. 



MUU « MU(JL» 



BD295 

500. 



SUM » 0. 



B0296 

501. 



Cl a 1, - bet^oy/muu 



BD297 

502. 



DO 1S4LJ * 1/8 



BD300 

503. 



PHASCO a PHASElLJ/JUl 



B0301 

504 . 



PHASKO a 2. » PHASCO 



BD302 

505. 


184 

SUM a SUM + A(LJ)*(PHASCO^fIP0S(LJ/J + 1> + 

PHaSKO+INEG(LJ/J + 

1 ) ) 

BD303 

506. 



TERM » INEg( JL/JJ 



BD304 

•brj*7 , 



TMFr < II . M- ■ Ti ' T|sl‘‘r? I II . 1 i i 1 

\ /Mini 


Rn*qn5 


uUo • 


IflU 

iA'-\All > i-A<vAii iiuUl UU/ w 

wU.^Uw 

509* 

C 




510. 

C 

TEST R0UTINE--CHECK FOR COMPLETION OR DIVERGENCE— PRINT IF REACHED 

' 

511. 

c 

convergence or on abort STATUS 


512. 

c 




513. 



IF( JK - 6)l91il91i650 

BD3061 

51*t. 


650 

IF(PARAM - 10.)191/190/190 

BD307 

515. 


191 

IF{PARAM - TOL ) I93i 192# 192 

BD308 

516. 


192 

IF( IK - 1 )l94il94#l95 

B0309 

517. 


195 

IF(AB0RT)2OO/2OOil96 

BD310 

518. 


190 

IF( AB0RT)71Oi7lO#196 

BD311 

519. 


710 

ABORT =1 

BD312 

520. 



PRINT 211#JK 

BD313 

521 . 


211 

FORMAT! lX///5Ki »¥¥4iPROQRAM HAS DIVERGED AFTER»#I3#» I TERAT I ONS— I N 

BD314 

522. 


IITIATE AB0RT<^^#5///) 

BD315 

523. 



GO TO 133 . i 

BD316 

524. 


193 

in IK - 1 )196#196#198 

B0317 

525. 


198 

PRTVAR » 10 

B0318 

526. 



GO TO 199 

BD319 

527. 


194 

IF!JK - 2)196/196/202 

BD320 

528. 


202 

IF(JK/PRINT - (JK + PRINT « 1 ) /PRINT 1 200/ 1 96/ 200 

BD321 

529. 


196 

PRTVAR B 1 ^ 

BD322 

530. 

c 




531. 

c 

PRINT IPOS . ... 


532. 

c 




533. 


199 

PRINT 500/JK 

BD323 

534. 



PRINT 80 . . 

60325 

535. 



DO 910JR » 1/NDS . ' 

00326 

536. 



Q ( JR ) = 0 . 

8D327 

537. 



DOR ( JR ) B 0* 

BD328 

538. 



DO 20SJQ B 2^6 

BD329 

539. 



JA = JQ 

BD3291 

540. 



IF(JR " LMTfJOl 1209/209/208 

BD330 

541 . 


2C8 

CONTINUE 

B0331 

542. 


209 

KFF B KF(KK/JA - 1) 

BD332 

543. 



DO 210JQ B 1/8 

B0334 

544. 



Q(JR) B (3(JR1 + A( JQ)4MU( J0)¥(1P9S(JQ/JR) - INEQ(JQ/JR)) 

BD335 

545. 


210 

DQR(JR) » DQR(JR) + A ( JQ ) POS ( JQ/ JR ) + INEQIJQ/JR)) 

BD336 

546. 



0( JR) = 0(JR)¥2. 43. 141593 

BD337 

547. 


910 

DaR(JR) B "2.43.1415934KFF4DQR( JR) + 4 . 4KFF4S4T ( JR ) 4444NN+F ( KK/ JR ) 

BD33S 

548. 



DO 207JR B 1/NDS/PRTVAR 

BD3381 

549. 



RJ B JR 

BD339 

550 « 



PRINT 281/RJ/ (IPOSI JS/ JR)/JS » 1/ 8 ) / 0 ( JR ) / DQR ( JR ) 

BD340 

551 . 


281 

FORMAT! lXiF4.0/2X/9(E10.5/ lX)/Eli.5) 

BD341 

552. 



IF!IK - 1)2200/2200/207 

BD342 

553. 


2200 

IFIPARAM - TeU2201/207/207 

603421 

554. 


2201 

IF((JR - D/IO - (JR + 81/10)207/2220/207 

BD343 

555. 


2220 

JT » JR/10 +1 

BD344 

556. 



SAVE! JT/1 I ■ JR 

0D345 

557. 



L « JT 

BD346 




56cr. 


SAVE( JTi 10) « Q( JR) 



BD349 

561 . 


SAVE( JT/ 11 ) » D0R( JR) 



B0350 

562* 


207 CaNTINUE 



B0351 

563* 


PRINT 82 



BD352 

sen- 

C 





ses* 

c 

PRINT INEG 




566* 

c 





567* 


Oe 215JR « NDSj li-PRTVAR 



BD353 

568* 


R J * J R 



BD354 

569. 


215 PRINT aSljRJM INEG( JSi JR)#JS • l/8j 



BD355 

570. 


0U7PUT PARAH 



BD3551 

571 * 


IFtJK - 1)200/200/629 



BD3552 

572* 


629 IF( IK - 1 1216/216/217 



B0356 

573. 


216 IFtPARAM - T0U 107/ 2 1 7/ 21 7 . . . 



BD3561 

574. 


107 D0 21SJR » NOS/ 1/-10 



BD357 

575. 


JT “ (NDS - JR + D/iO + 1 + t 



BD358 

576. 


SAVE! JT/1 ) = JR 



BD359 

577. 


SAVE! JT/10) = 0. 



BD360 

578. 


SAVE! JT/11 ) “ 0. 



BD361 

579. 


00 218JZ - 2/9 



B0362 

580. 


218 SAVEiJT/JZ) = INEQIJZ - 1/JR) 



BD363 

581 • 

c 





582. 

c 

.CHECK F0R NEED 0F N0ND1MENSIONAUIZINQ-..SKIP 

SECTION IF 

ON abort 

STATUS 

583. 

c 





584. 


217 IF(PARAM - 10.1220/224/224 



BD364 

585. 


224 IF(JK - 6)220/220/2000 



BD3641 

586. 


220 IF( ABORT I221/221i200 



BD365 

587. 


221 IF{N0NDM1222j222/223 



BD366 

588. 


223 1F(PARAM - TOL) 399/ 200/ 200 



BD3661 

589. 

c 





590. 

c 

dimensionless RESULTS--SAVE CALCULATIONS IF 

ON SPECIAL 

band 


591. 

c 





592. 


399 Cl « l./IPOSN 



BD367 

593. 


C2 » 1./Q00 



BD368 

594. 


C3 « l./KF(KK/5)/Q0O 



BD369 

595. 


PRINT 86 



. BD370 

596. 


PRINT 80 



BD3701 

597. 


DO 225JR =» 1/NDS/PRTVAR . 


! 

B0371 

598. 


DMLQR 0(JR)*C2 


J 

BD372 

599. 


DMLDQR » DQRtJRl»C3 



BD373 

600. 


DO 226JP » 1/8 



BD374 

601 . 


226 OML(JP) = IPOS! JP/JR)¥C1 



BD375 

602* 


IF( IK - 1 )227/2a7/228 



BD376 

603. 


227 IF(PARAM - TOL ) 1 09/ 228/ 228 



3761 

604. 


109 IF((JR - 1)/10 - UR ^ 8 )/10)228/229/228 



377 

605 • 


229 JP » JR/10 + 1 + L + L 



BD378 

606 • 


SAVEUP/1 ) » JR 



6D379 

607. 


SAVE ( JP/ 10 ) DMLQR 



0D38O 

608* 


SAVEt JP/ 11 ) = DMLDQR 



BD381 

' n - 





nn'’ « 


i-l. 


w i U • 


U w)J 


W u V.' 

611% 


228 

RJ * JR 

B0384 

612* 


225 

PRINT 281>RJ^ DMt/DMLQR/DMLDQR 

BD385 

613. 



PRINT 82 

BD386 

614* 



00 254JR = NDSi 1#-PRTVAR 

BD387 

615. 



RJ = JR 

. . . 60388 

616. 



D0 251JP « 1#8 

BD389 

617. 


251 

DMUJP) « INEGl JP/ JR)»C1 

B0390 

618. 



IFtlK - 1)252/252/254 

BD391 

619. 


252 

IFtPARAM - T0U57/254/254 

BD391 1 

620. 


57 

IF((JR - D/IO « <JR + 8^/10)254/253/254 

BD392 

621 . 


253 

JP - INDS - JR + D/IO + 1 + U 4 L + t 

BD393 

622. 



SAV£{ JP/1) - JR 

BD394 

623. 



SAVE(JP/10) “0. 

BD395 

624. 



SAVE(JP/11) “0. 

BD396 

625. 



D0 255JT - 2^9 

1 BD397 

626 . 


255 

SAVE! JP/ JT) « DML( JT - 1 1 

BD398 

627. 


254 

PRINT 281/RJiDMU 

BD399 

628. 


222 

IF(PARAM ^ T0U26O/2OO/2OO 

BD400 

629. 


200 

CONTINUE 

BD401 

630. 

C 




631. 

C 

IF 

PROGRAM PROGRESSES TO THIS POINT/ IT TO0K 

TOO many ITERATIONS--ABORT 

632. 

C 




633- 



IF( AB0RT)27O/27Oi2OOO 

BD402 

634. 


270 

ABORT = 1 

BD403 

635. 



PRINT 271/KK 

BD404 

636. 


271 

FORMATUX///’ 4>f4BAN0»/ 13/ » WENT BEyONO 

MAXIMUM ITERAT IONS"- INI TI BD405 

637. 


1AT2 ABPRT4^¥»///) 

BD406 

638. 



G0 T0 .133 

, , BD407 

639. 

c 




640. 

c 

STARTING ROUTINE F0R INEGP# 0R REFtECTANCE CAUCULATI0NS 

641 . 

c 




642. 


260 

IF ( NN - 1 ) 732/732/ 733 

BD4071 

643. 


732 

00 734JP - 1/8 

BD4072 

644. 


734 

INEGP(JP/U » IN£G(JP/1) 

B04073 

645. 



GO TO 360 

BD4074 

646 . 


733 

PRINT 351 

BD4075 

647. 



DO 285JP » liNOS 

BD408 

64S. 



J « 0 

■ 1 BD409 

649. 



A 1 » 0 * 

BD410 

650. 



DO 280KJ ..8/1/*! 

BD411 

651. 



IF(MU(KJ1 " 017)283/283/282 

BD412 

652. 


282 

J « J + 1 

BD413 

653. 


280 

A1 « A1 + INEQ(KJ/JP) 

BD414 

654. 


283 

in J)2810/2810/.284 

B0415 

655. 


234 

Al » Al/J 

BD416 

656 t 



GO TO 58 

BD417 

657. 

2810 

Al - INEQ(S/JP) 

BD4171 

658. 


58 

00 285KJ « 1/8 

BD4172 

659. 


285 

INEGP(KJiJP) « Al 

B0418 

r» 



-^nn i*-' i . 

n 1 « 



662V 

C 



663* 

c 

FRONT WALL CALCULATIONS 


664* 

c 

- 


665* 


D0 3011 » 1/S 

BD421 

666 • 
667 • 

c 

301 IPeSPdill » RFL8( n4lNEGP( I/l ) f (1* ^ RFLO ( I) ) ♦ IPOSN^NN ^ 

B0422 

668* 

c 

MARCH T0 back wall 


669. 

c 



670. 


D0 3101 =1/5 

BD423 

671. 


D0 3031L » 1/8 

BD424 

672* ' 


D0 303IN » 1/8 

B0425 

673. 


303 PHASE(IL/IN) » 1* + I SOT (1 J »MUP ( IL 1 4MU ( IN ) 

BD426 

674. 


III « LMT(I » + 1 

B0427 

675. 


JJJ. = LMT( I + 1 ) 

BD428 

676 . 


IF( III - 1)304,304,305 ) 

B0429 

677. 


304 III - 2 

BD430 

678- 


305 SFF « SFIKK/ I) 

80431 

679. 


KFF = KF ( KK/ I ) 

BD432 

680* 


SET « SFF + KFF 

BD433 

681 . 


C2 * DY^SFF/2. ... 

BD4331 

682. 


C4 « DY^S-^KFF/3. 1415934NN 

BD4332 

683. 


D0 310J. « III/JJJ 

BD434 

684. 


C3 ■ C4^r(J “ 1}¥#4+F(KK/J) 

B04341 

685. 


00 310JL "1/8 

B0435 

686 . 


MUU = MUPlJU 

BD436 

687. 


SUM « 0. 

BD437 

688. 


Cl = 1. - BET4DY/MUU 

BD438 

689. 


00 3ULJ » 1,8 

BD441 

690. 


PHASCe * PHASEUL/LJl ... 

BD442 

691 . 


PHASK0 « 2. • PHASCO 

BD443 

692. 


311 SUM “ SUM + A (LJ)4(PHASCO4lP0S(LJ#J-il + PHASK04 I NEQ t L J/ J-1 ) ) 

BD444 

693. 


TERM - IP0SP(JL/J) 

BD445 

694. 


IPOSPIJL/J) « C14IP0SP( JL/ J - 1) + IC24SUM + C3I/MUU 

BD446 

695. 


310 PARAM = PARAM + (TERM/IP0SN - I P0SP 1 JL/ J ) /IPOSN ) 

' B0447 

696. 

c 



697. 

698* 

c 

c 

BACK eaUNDARY CONDITIONS 


699. 


IF(KFRES - 1 )118/118/119 . , . : 

BD4471 

700. 


119 Cl • S^NN*T(NDS/4 + 4/3.141593i^FIKK/J) 

004472 

701 . 


00 311 » 1/8 

BD4473 

702. 


31 INESP(I/N0S) = FRFLOn )4lP0SPd/ND5) II. - RBBJ4CI 

B04474 

703. 


G0 T0 800 

BD4475 

704. 


118 Cl « (1. - RBB)4S^NN*T(N0Sl444/3»l'fl593»F(KK/J) 

BD443 

705. 


00 3121 = 1/8 . , 

. B0449 

706. 

707. 

c 

312 INE^Pd/NDS) » R8B»IP0SPd/NOS) .4 Cl 

BD450 

70S. 

c 

MARCH T0 FRONT WALL 


709. 

c 



710. 


800 D0 3151 » 5/1/-1 

BD451 


14 - 

15 « 

16 ^ 

\7< 

'18 

'19 

'20 

'21 

'22 

'23 

'24 

'25 

'26 

'27 

'28 

'29 

'30 

■31 

'32 

'33 

'34 

'35 

'36 

'37 

■38 

'39 

'40 

■41 

'42 

'43 

'44 

'45 

'46 

'47 

'48 

'49 

■50 

51 

52 

53 

54 

55 
'56 
57 
'58 
'59 
60 
61 


C 

c 

c 


c 

c 

c 


UU-wAwJ. iJU ■ ^ ^ 

316 PHASE(IL/IN) * 1. * IS8T ( I ) #MUP * lU 4MU ( IN ) BD454 

III « LMT( I + 1 } BD455 

JJJ » UMT( I ) + 1 BD456 

IFdll - NDSl3l7^3l8i318 BD457 

318 III NOS *1 B0458 

317 SFF « SF(KK/I) BD459 

KFF ■» KF ( KKj I ) BD460 

SET « SFF + KFF B0461 

C2 * DY*SFF/2. - BD4611 

C4 a DY¥S¥KFF/3.14l593*NN BD4612 

00 315J a III/JJJ/-1 ' BD462 

C3 » C4»^T(J + i l**4»fF(KKi J» BD4621 

00 315JL a 1^8 BD463 

MUU a MUP(JL) BD464 

SUN a 0* I B0465 

Cl a !• - BET¥DY/NUU BD466 

DO 320LJ a 1,8 BD469 

PHASCe a PHASE(JLiLJ) BD470. 

PHASK6 a 2. - PHASCO BD471 

320 SUN a SUN + A (LJ)'MPHASCe*IP0S(.LJjJ+l) + PHASK0* I NEQ (L J/ J + 1 ) I BD472 

TERM a INEGPULiJI 60473 

INEGP(JU/J) a C1 *INEGP(UUj J + 1* + (C2*SUM + C3»/NUU BD474 

315 PARAM a PARAN + tTERM/IPOSN - INEGP ( JL/ U 1 / IP0SN ) 442 BD475 

CHECK F0R divergence^ S0UUTI8N, 6R AB6RT STATUS— -ST ANOARD PROCEDURE AS ABOVE 

1F(JP - 1 )323i323i324 BD476 

324 IFIPARAN a 1 0 I 323 j 325, 325 BD477 

323 IFtPARAN - TOL) 794, 327, 327 BD478 

325 IF(AB0RT)33O, 330, 2000 ' BD479 

327 1F(AB0RT)34O,34O,321 B04791 

330 ABORT =1 BD480 

PRINT 331, JP BD481 

331 FORMATdX///' 44¥REFLECTANCE CALCULATION DIVERqED--INITI ATE ABORT BD482 

1 AFTER', 13,' ITERATI0NS444 V//* B0483 

GO TO 260 B0484 

340 IF(IK - 1)321,321,300 B0485 

794 IFUK - 1)321,321,326 . , BD4851 

326 PRTVAR a 10 ! 00486 

GO T0 350 60487 

321 PRTVAR * 1 BD48R 

\ 

PRINT REFLECTANCE FIELD 

350 PRINT 500, jP B0489 

351 FORMAT; IX///' REFLECTANCE FIELD'//) . ■ BD490 

PRINT 352 BD491 

352 F0RNAT(' NODE ' , 6X, ' 1 1 ' , 9X, ' 1 2 ' , 9X, ' 1 3 * , 9X# * 1 4 t , 9X , ' 1 5 ' , 9X, ' 1 6 ' , 9X, BD492 

1 * 17 ' ,9X, ' 18 ' // I B0493 

nfl “‘t=i'5T - 1 . • '^'■'TV * R Rr»-q' 


/ 


/ w 


IV J, i • UU T ^ J 


76 41 


353 

PRINT 281/RI# 1 1P0SP( IX# I IX ■* 1#8) 

BD496 

765. 



PRINT 82 

604961 

766. 



DO 3541 » NDSil#-PRTVAR 

B04962 

767. 



RI - I 

BD4963 

768. 


354 

PRINT 281#RIi ( INEGPUX# I »/IX « IjS) 

BD4964 

769. 



OUTPUT PARAM 

BD4965 

770. 



IFIPARAM - T0L)36O#3OO#3OO 

B0497 

771 . 


300 

CONTINUE 

80498 

773. 



IF ( ABORT l370#370/2000 

BD499 

773. 


370 

ABORT X 1 

BD500 

774. 



PRINT 371 , 

80501 

775. 


371 

FORMATS IX///' *^»REFLECTANCE ITERATIONS EXCEEDED MAXIMUM—INITIAT 

BD503 

776. 



IE AB0RT¥»4'///) 

BD504 

777. 



GO TO 260 . 

BD505 

778. 

C 


i 


779. 

C 

CALCULATE REFLECTANCE FOR BAND AND PRINT 


780. 

C 




781 . 


360 

R(KK) =0 

BD506 

782. 



DO 3611 = i#e 

BD507 

783. 



RF(I) » RFLO^n + iNEGPt I/U/NNJ»ta, - RFL0 ( I ) ) /IPOSN 

BD50S 

784. 


361 

R(KK) X R(KKJ + AUl ( I >+NUl ( I )4HRF( I) 

BD509 

785. 



PRINT 921#RF 

BD5091 

786. 


921 

FORMATCIX//' REFLECTANCE FIELD! »//7X#B(F10.8#lX)/) . 

BD5092 

787. 



R(KKI * R(KK)/. 50151552 

0D5O93 

788. 



OUTPUT R(KK) 

BD510 

789. 

C 




790. 

c 

SAVE QR AND DORDY FROM THIS BAND FOR OVERALU RESULTS 


791 . 

c 




792. 



DO 3621 = 1#NDS 

BD511 

793. 



QRSAVl I ) B QRSaV( I } + Qt I 1 

BD512 

794. 


362 

OQRSAVU) » DORSAVtl) + DORdl 

8D513 

795. 



OARFLl « OARFLl + IPOSN^RtKK) , , 

BD514 

796. 



0ARFL2 a 0ARFL2 + IPOSN 

B0515 

797. 


65 

CONTINUE 

BD516 

798. 

c 




799. 

c 




800* 

c 

CALCULATE OVERALL REFLECTANCE AND PRINT# ALONG WITH OVERALL OR AND DORDY 

801 . 

c 




802. 



0ARFL =. OARFL1/0ARFL2 ' 

BD517 

803. 



PRINT 410iOARFL# (QRSAVI Dil « 1#NDSI 

BD518 

804. 


410 

F9RNAT('1 over all reflectance -*#F8.7///* OVERALL FLUXES# STARTI 

BD519 

805. 



ING AT FRONT W ALL ! ' // ( 1 X# lOE 1 2 . 5/ H 

BD520 

806. 



PRINT 411# { DQRSAV 1 n# I ■ liNDSI 

B0521 

807. 


411 

FORMAT! IX/// ' OVERALL D I VERGENCE! » // ( IX# itfE 1 2 • 5/ ) ) 

8D522 

808. 

c 




809* 

2000 

END 

BD523 





Test C,asb 1. 

to= 1.0 

Rg=.os 

X - o Qrsorf^Ptc. Scat7-e:«.//jg- 3 

J^fOT>€X “ 1'^ 

l.z FO^ Z^ZAaJT^ 



NOS = 101 

TQL « 9099953E-11 
TEST » 1000 
PRINT = 25 
N8NDM » 1 

THICK « 1*00000 CM 
BNDS » IIOOOOOOOO 

SPECIAL -- 0 STANDARD -- 1 2 


TEMPO » 1 

TEMP « 

*000 

*000 

• 000 

• 000 

• 000 

• 000 

• 000 

• 000 

• ooo 

• ooo 

(KELVIN) 


• 000 

• 000 

• 000 

• 000 

• 000 

• 000 

• 000 

• 000 

• ooo 

• ooo 


.000 

• 000 

• 000 

• 000 

• 000 

• 000 

.000 

• 000 

• ooo 

• ooo 


• 000 

• 000 

• 000 

• 000 

• 000 

• 000 

• 000 

' •ooo 

.000 

• ooo 


• 000 

.000 

• 000 

• 000 

• 000 

*000 

• 000 

• 000 

• 000 

• ooo 


• 000 

• ooo 

• 000 

• 000 

.000 

• 000 

.000 

• ooo 

• ooo 

.000 


• 000 

.000 

• 000 

• 000 

.000 

• 000 

• 000 

• ooo 

• ooo 

• ooo 


• 000 

• 000 

• 000 

• 000 

• 000 

• 000 

• 000 

• ooo 

.000 

• ooo 


• 000 

.000 

• 000 

• 000 

• 000 

*000 

• 000 

• ooo 

• ooo 

.000 

IS0T - 

• 000 
• 000 

•0000000 

.000 

•0000000 

•000 *000 
•0000000 • 

• 000 
0000000 

.000 

•0000000 

• 000 

.000 

.000 

.000 


INTENSITIES IN WATTS/CM^¥2/STERADI AN 

FLUXES IN WATTS/CM¥42 

FLUX DIVERGENCE IN WATTS/CrU^^^S 



BAND 1 


INDEX • l.AOO 
RB » t030 

Q0 » 100.0000000 WATTS/CN^^2 

K « <000000 .000000 <000000 <000000 <000000 

S »» 1<000000 1<000000 l<000000 1<000000 ItOOOOOO 


GAUSSIAN INTEGRAL 


MU»S 

<16699863 

•4627501A 

<65259266 

<71393394 


WEIGHTS 

<32747149 

<25248051 

<11990213 

<03555630 


<76911736 

.84992695 

<93073654 

<98591995 


<07182908 

<08537477 

•07182908 

<03555630 



WAVELENGTH INTERVALS <OOOOOE 00 - .99999E 05 CM 

<OOOOOE 00 - <99999E 09 MICRONS 


F(0 - LT) ; 


• 76529E 

-80 

<4nS0E-82 

. 10295E-83 

*51476E-84 

•OOOOOE 00 

•OOOOOE 

00 

•OOOOOE 00 

•OOOOOE on <00000E 00 

•OOOOOE 

• OOOOOE 

00 

•OOOOOE 00 

•OOOOOE 00 

•OOOOOE 00 

•OOOOOE 00 

•OOOOOE 

00 

•OOOOOE 00 

•OOOOOE 00 <OOOOOE 00 

• 0000(| 

.OOOOOE 

00 

•OOOOOE 00 

•OOOOOE 00 

•OOOOOE 00 

•OOOOOE 00 

•OOOOOE 

00 

•OOOOOE 00 

•OOOOOE 00 <5l476E-84 

• 6588^- 

.32680E 

-79 

•26313E-77 

•62500E-01 

•OOOOOE 00* 

•65446E-77 

•OOOOOE 

00 

•80964E-77 

.22204E-15 <13092E-79 

• 5037W* 

•OOOOOE 

00 

•OOOOOE 00 

•43959E-72 

.86736E-18 

<22041E-38 

•46137E 

08 

•35112E 51 

i 

•13605E 46 .a5759E 51 

• 1880^- 

• 22859E 

09 

.40U6E-17 

.27636E-75 

•27537E 48 

.44060E 49 

•94040E-37 

•42409E-25- 

•21828E-09 <25760E 51 

•18877E- 

•27538E 

48 

•23510E-37- 

.96207E 12 

•78603E-21 

•27638E-75 

•44061E 

49 

•13346E 46 

•60818E 08 <18809E-36 

• 6112^ - 

.35112E 

51 

•43656E-10 

. 11755E-37 

•36734E-39 

•12148E-26 

•OOOOOE 

00 

•25593E-54 

•48172E-37 .10408E-16 

*29387^- 

.33144E 

13 

•35903E-73 

.36734E-39 

• 18808E-36- 

.21690E-06- 

•24501E-05 

•19966E 46 

•43369E-17 .44Q65E 49 

• 1410«- 

<60S1SE 

08 

.25760E 51 

.27541E 48 

•94039E-37 

•20881E-52 

.25759E 

51 

•55039E-19 

•13306E 46-.56548E*10 

-•9244a- 



/ ^ J w W w 


:wrOG « 45.584V 

RFLI 1.0000000 1*0000000 

1.0000000 

.3219466 

•Q909002 

RFL0 • .5718594 *3050691 

•0863307 

•0465407 

•0332410 


ITERATION 19 









NODE 

11 


12 


13 


14 


15 

1. 

• 23624E 

02 

•18984E 

02 

•16367E 

02 

•47356E 

02 

.58092E 

11* 

• 24261E 

02 

.20172E 

02 

. 17607E 

02 

.44420E 

02 

•54038E 

^1 • 

.24142E 

02 

•20930E 

02 

. 18528E 

02 

.41740E 

02 

•50359E 

31 . 

• 23537E 

02 

.213UE 

02 

. 19150E 

02 

•39258E 

02 

•46988E 

'4l. 

.22608E 

02 

.21361E 

02 

•19497E 

02 

•36931E 

02 

•43870E 

01. 

.21455E 

02 

.21124E 

02 

. 19592E 

02 

•34725E 

02 

•40962E 

61. 

.20136E 

02 

•20637E 

02 

.19457E 

02 

•32610E 

02 

•38225E 

/I • 

. 18681E 

02 

. 19929E 

02 

.19111E 

02 

•30562E 

02 

•35627E 

31. 

. 17099E 

02 

.190H2E 

02 

• 1856SE 

02 

•28556E 

02 

•33136E 

31 . 

•15380E 

02 

. 17925E 

02 

•17834E 

02 

•26568E 

02 

•30720E 

101 • 

•13489E 

02 

•16634E 

02 

•16908E 

02 

•24567E 

02 

•28346E 


101. 

•40467E 

00 

• 49901E 

00 

•50723E 

00 

•73702E 

00 

•85037E 

91. 

•54139E 

01 

.26000E 

01 

.20362E 

01 

•21123E 

01 

•21185E 

31. 

•90795E 

01 

• 470UE 

01 

.36473E 

01 

•35833E 

01 

.34S91E 

7i . 

•11922E 

02 

.67598E 

01 

•52977E 

01 

•510S9E 

01 

•49230E 

61. 

•14252E 

02 

•87547E 

01 

•69602E 

01 

•66621E 

01 

•63938E 

61* 

•16252E 

02 

• 10676E 

02 

.86163E 

01 

.S2238E 

01 

•78827E 

41. 

•18028E 

02 

.12517E 

0? 

•10253E 

02 

.97797E 

01 

•93751E 

31 . 

.1963SE 

02 

. 14276E 

02 

.11858E 

02 

.11318E 

02 

•10859E 

21. 

.21109E 

02 

•15947E 

02 

•13421E 

02 

.12S26E 

02 

•12321E 

u. 

•22443E 

02 

•17521E 

02 

•14930E 

02 

•14291E 

02 

•13749E 

1. 

PARAM 

•23624E 02 
« 8.441262E 

.1B984E 

-12 

02 

.16367E 

02 

•15696E 

02 

•15126E 


DIMENSIONLESS 









NODE 

11 


12 


13 


14 


15 

1* 

• 74218E 

00 

.59641E 

00 

.51417E 

00 

.14877E 

01 

•18250E 

u. 

• 76217E 

00 

•63371E 

00 

•55316E 

00 

•13955E 

01 

•16977E 

. 

. 7'^R 

on 


on 


00 

• 1 31 

0 J 

. 1582 IP 


•0394923 

.0292094 

•0278213 

.0290592 

.0279618 

•0277839 








l)dto>y 

o. 



16 


17 


18 


OR 



02 

.60491E 

02 

•60965E 

02 

•61020E 

02 

•71419E 

02 

^0 

02 

•56534E 

02 

•57286E 

02 

•57533E 

02 

•71445E 

02 

-■ 0 

02 

•52908E 

02 

•53883E 

02 

•54288E 

02 

•71466E 

02 

^ 0 

02 

•49556E 

02 

•50708E 

02 

•51243E 

02 

•71482E 

02 

» 0 

02 

.46433E 

02 

•47725E 

02 

•48368E 

02 

•71495E 

02 

■ 0 

02 

•43501E 

02 

•44903E 

02 

•45635E 

02 

•71507E 

02 

■ 0 

02 

.40728E 

02 

•42216E 

02 

•43020E 

02 

•71519E 

02 

. 0 

02 

•38084E 

02 

•39639E 

02 

•40503E 

02 

•71532E 

02 

aO 

02 

•35543E 

02 

•37149E 

02 

•38062E 

02 

.71547E 

02 

■ 0 

02 

•33077E 

02 

•34722E 

02 

•35676E 

02 

•71568E 

02 

0 

02 

.30653E 

02 

•32331E 

02 

•33320P 

02 

•71597E 

02 


00 

•91960E 

00 

•96992E 

00 

•99961E 

00 




01 

•20655E 

01 

•20158E 

01 

•19866E 

01 




01 

•33184E 

01 

•31700E 

01 

•30819E 

01 




01 

•46419E 

01 

•43992E 

01 

•42538E 

01 




01 

.60114E 

01 

•56802E 

01 

• 54804E 

01 




01 

•74087E 

01 

•69957E 

01 

•67449E 

01 




01 

•88194E 

01 

•83320E 

01 

•80340E 

01 




02 

.10231E 

02 

•96768E 

01 

•93358E 

01 




02 

.11632E 

02 

•11018E 

02 

•10639E 

02 




02 

• 13008E 

02 

•12343E 

02 

.11930E 

02 




02 

• 14343E 

02 

•13636E 

02 

*131932 

02 






16 


17 


18 


OR 



01 

.19004E 

01 

. 19153E 

01 

•19170E 

01 

•71419E 

00 

■o 

01 

• 17761E 

01 

•17997E 

01 

•18075E 

01 

•71445E 

00 

■ r 

01 

. 1 6*??E 

01 

•16928F 

01 

. 17055P 

01 

•71466E 

00 

• c 


w i. • 

• / ^ J u- 

uu 

. U W vJ U t.. 

VJ.U 

. ‘JU 1 UW1_ 

uu 

' i C- ^ v> t_ 

Kt 1 

•It/ W pw 1_ 

U 1 

• A w ^ w 

U A 

« A <J ^ V A (-■ 

<J A 

• A ^ ^ . 

U A 

• . * c- 

SJ . 

41. 

•7i026£ 

00 

•67107E 

00 

•61251E 

00 

. 11602E 

01 

•13782E 

01 

.145S7E 

01 

• 14993E 

01 

. 15195E 

01 

•71495E 

00 

5^1 . 

«^67403E 

00 

•66362E 

00 

.61549E 

00 

. 10903E 

01 

. 1E869E 

01 

. 13666E 

01 

. 14107E 

01 

•14337E 

01 

•71507E 

00 

61. 

•63259E 

00 

•64832E 

00 

.61126E 

00 

. 10245E 

01 

.12009E 

01 

•12795E 

01 

.13263E 

01 

•13515E 

01 

•71519E 

00 

71. 

•58687E 

00 

•62609E 

00 

•60039E 

00 

•96014E 

00 

•11193E 

01 

.11964E 

01 

•12453E 

01 

.12724E 

01 

•71S32E 

00 

81. 

.53718E 

00 

•59758E 

00 

.58333E 

00 

•89713E 

00 

•10410E 

01 

. 11166E 

01 

. U 6 7 1 E 

01 

. 11958E 

01 

.71547E 

00 

91. 

.4S318E 

00 

•56312E 

00 

.56028E 

00 

•83466E 

00 

.96510E 

00 

•10391E 

01 

.10908E 

01 

.1120BE 

01 

.71568E 

00 

101 . 

•42377E 

00 

.52256E 

00 

.53117E 

00 

•77181E 

00 

•89050E 

00 

•96300E 

00 

•10157E 

01 

•10468E 

01 

•71597E 

00 


iOl. 

•12713E- 

01 

.15677E 

-01 

. 15935E-01 

.23154E 

•01 

•26715E-01 

.28890E-01 

.30471E- 

01 

•31404E 

■01 

91. 

. 17008E 

00 

.81681E 

-01 

.63970E-01 

•66359E 

-01 

•66556E-01 

•64891E-01 

.63328E- 

01 

.6241 IE 

-01 

81. 

.28524E 

00 

. 14769E 

00 

. 11458E 

00 

• U257E 

00 

•10961E 

00 

. 10425E 

00 

•9958SE- 

01 

.96819E 

■ 01 

71. 

•37455E 

00 

.21236E 

00 

. 16643E 

00 

.16050E 

00 

.15466E 

00 

•14583E 

00 

. 13820E 

00 

•13364E 

00 

61. 

. 44774E 

00 

•27504E 

00 

.21866E 

00 

•20930E 

00 

.20087E 

00 

.18885E 

00 

•17845E 

00 

•17217E 

00 

51. 

.51058E 

00 

.33538E 

00 

•27069E 

00 

.25836E 

00 

•24764E 

00 

•23275E 

00 

.21978E 

00 

.21190E 

00 

41. 

.56637E 

00 

.39324E 

00 

.32209E 

00 

•30724E 

00 

•29453E 

00 

.27707E 

00 

.26176E 

00 

.25240E 

00 

31. 

•61695E 

00 

.44S49E 

00 

.37253E 

00 

.35555E 

00 

.34113E 

00 

.32142E 

00 

•30401E 

00 

.29329E 

00 

21. 

.66315E 

00 

.50098E 

00 

•42164E 

00 

.40294E 

00 

.38708E 

00 

.36542E 

00 

.34615E 

00 

.33422E 

00 

11. 

• 70506E 

00 

.55044E 

00 

•46903E 

00 

. 44S96E 

00 

•43194E 

00 

. 40865E 

00 

.38778E 

00 

.37478E 

00 

1. 

.74218E 

00 

.59641E 

00 

.51417E 

00 

.49311E 

00 

•47518E 

00 

.45060E 

00 

.42839E 

00 

•41448E 

00 

reflectance field 
















ITERATION 4 








sOOE 11 

12 

13 

14 

15 

16 

17 

18 


1. 

•35743E 

02 

.52783E 

02 

.58305E 

02 

•60162E 

02 

•60781E 

02 

.609712 

02 

•61016E 

02 

.61020E 

02 

11. 

. 34315E 

02 

. 49200E 

02 

•54242E 

02 

•56141E 

02 

•56937E 

02 

•57304E 

02 

•57482E 

02 

•57562E 

02 

21. 

.32944E 

02 

•45952E 

02 

-50553E 

02 

.52466E 

02 

•53399E 

02 

•53910E 

02 

.5419SE 

02 

.54341E 

02 

31. 

•31601E 

02 

.42972E 

02 

.471722 

02 

•49079E 

02 

•50116E 

02 

•50743E 

02 

•51121E 

02 

.51317E 

02 

41. 

.30262E 

02 

. 40208E 

02 

. 44045E 

02 

•459302 

02 

•47046E 

02 

•47766E 

02 

.48219E 

02 

•48459E 

02 

51'. 

.28914E 

02 

.37619E 

02 

.4U28E 

02 

•42980E 

02 

•44155E 

02 

.44949E 

02 

.45463E 

02 

•45739E 

02 

61. 

•27545E 

02 

.35170E 

02 

•38382E 

02 

.40194E 

02 

•41412E 

02 

•42266E 

02 

.42830E 

02 

•43136E 

02 

71. 

.26143E 

02 

•32829E 

02 

.35776E 

02 

.37544E 

02 

•38791E 

02 

•39691E 

02 

.40297E 

02 

.40629E 

02 

81. 

.2469SE 

02 

•30567E 

02 

•33277E 

02 

•34999E 

02 

•36266E 

02 

•37204E 

02 

•37843E 

02 

•38197E 

02 

91. 

• 23196E 

02 

.28356E 

02 

•30854E 

02 

.32531E 

02 

•33811E 

02 

•34779E 

02 

.35447E 

02 

•3581RE 

02 

101 • 

.21617E 

02 

.26161E 

02 

.28473E 

02 

•30108E 

02 

•31397E 

02 

.32389E 

02 

.33082E 

02 

•33468E 

02 


ul . 

» 64o30L 

00 

• /u48dt 

uu 

.ot>419t 

00 

. 9UJ24E 

00 

*^41ylt 

00 

.9/ 168t 

00 

• 99c4t>t 

00 

• 10040E 

01 

91. 

•20554E 

01 

.21286E 

01 

•21166E 

01 

•20806E 

01 

.20438E 

01 

.20140E 

01 

• 19935E 

01 

• 19824E 

01 

81. 

'•35552E 

01 

•35704E 

01 

.34817E 

01 

.33641E 

01 

.32534E 

01 

.31647E 

01 

•31030E 

01 

•30689E 

01 

71. 

. 5 1 0 66E 

01 

•50696E 

01 

.49105E 

01 

.47166E 

01 

•45358E 

01 

•43905E 

01 

•42888E 

01 

•42323E 

01 

61. 

•66824E 

01 

•65992E 

01 

•63766E 

01 

•61129E 

01 

•58670e 

01 

•566S3E 

01 

.55285E 

01 

.54507E 

01 

51. 

•82639E 

01 

•81403E 

01 

•78614E 

01 

.75346E 

01 

.7229GE 

01 

•69809E 

01 

•68055E 

01 

•67075E 

01 

41. 

•98369E 

01 

•96783E 

01 

.93502E 

01 

•89673E 

01 

•86078E 

01 

•83144E 

01 

.81061E 

01 

•79894E 

01 

31* 

•11389E 

02 

. 11201E 

02 

• 10831E 

02 

• 10399E 

02 

.99910E 

01 

•96567E 

01 

•94185E 

01 

•92846E 

01 

21. 

• 12910E 

02 

. 12697E 

02 

•12290E 

02 

•11816E 

02 

•11366E 

02 

.10996E 

02 

. 10731E 

02 

•10581E 

02 

11. 

. 14385E 

02 

•14151E 

02 

•13716E 

02 

•13206E 

02 

•12721E 

02 

•12319E 

02 

•12030E 

02 

.11867E 

02 

1* 

ARAM 

.15798E 02 
» .000000 

•15549E 

02 

.15091E 

02 

•14553E 

02 

•14039E 

02 

•13610E 

02 

•13301E 

02 

.13126E 

02 


reflectance FIELD! 

i 

•68030000 .40318537 *30733223 .26895428 .25078058 .24087214 .23519987 -23233533 

R(KK| « .289466 




} 


BAND 2 


INDEX « i*200 

o *030 

ae - 100*0000000 wATrs/CM^#2 

; » *000000 .000000 *000000 

3 ■ 1*000000 1.000000 1.000000 


•000000 

1.000000 


toooooo 

1.000000 CH »#-1 


3AUSSIAN 1 
MU'S 

.V.13190150 
^.36549670 
,0.51544112 
-.0*57374948 
.^65597504 
J.77638483 
^*89679462 
«. 97902018 


NTEGRAL 

WEIGHTS 

*25864875 

* 19941813 
.09470302 
*05293041 
*10702837 
*12721205 

* 10702837 
*05298041 


•iAVELENGTH INTERVAL; *99999E 05 

.99999E 09 


.99999E 70 CM 
.99999E 74 MICRONS 


F(0 " LT) : 

.80289E-77~*23133E«06 

.79437E-80 

•OOOOOE 

00 

•OOOOOE 00 

• 41454E 

-75 

•OOOOOE 

00 

•OOOOOE 00 .OOOOOE 00 

•OOOOOE 

1 

.OOOOOE 00 *00000E 00 

•OOOOOE 00 

•OOOOOE 

00 

•OOOOOE 00 

.OOOOOE 

00 

•OOOOOE 

00 

•OOOOOE 00 .OOOOOE 00 

•ooooo| 

•OOOOOE 00 .OOOOOE 00 

•OOOOOE 00 

•OOOOOE 

00 

•OOOOOE 00 

•OOOOOE 

00 

•OOOOOE 

00 

•OOOOOE 00 tioEgsE-sa 

•13126E 

•33603E-79 .86360E-77 

.72370E 76 

•OOOOOE 

00- 

•67470E-79 

•OOOOOE 

00 

.20916E- 

77 

•44409E-15 .13046E-79- 

•44645! 

•OOOOOE 00 .OOOOOE 00- 

.44645E-07 

• 19259E- 

33 

•10486E OS 

•S5S55E 

43 

•16521E 

1 

46 

•13331E 46 *25760E 51 

•94750I 

•18011E 46 .70530E-36 

* 13Sa*tE-26 

.47742E- 

37 

.47742E-37- 

•96207E 

12 

.11102E- 

13 

•56424E-36 .4774EE-37 

M 

•25759E 

1 

•7346SE-39 .27538E 48- 

•77151E-65 

•58676E 

49 

•178aSE-65 

• 2.a769E 

■17 

•55022E-19 

•13532E 46 .65828E-36 

.35265] 

.25771E 51 .104082-16 

.58677E 49 

•35900E- 

73 

•13389E 46 

.80297E 

1 

00 

0 

1 

•OOOOOE 

00 

•94040E-37 .22041E-38 

•44064E 

•5777SE-33 *35900E-73 

* 10S42E-18 

•13346E 

46- 

*83266E-26 

•13676E 

46 

• 19966E 

46 

•60838E-17 .94039E-37 

•188791 

•25759E 51 .47091E-36 

* 18563E-17 

•41361E- 

55 

.47742E-37 

.47740E 

-37 

•27536E 

48 

•24686E-48- .24536E-05 

.302U| 




.f\ f • 

:RTDG = 
■FUI 

W j A 1 J. 

56.4437 

1.0000000 

1.0000000 

1.0000000 

.2648743 

.0530140 

.0153192 

*0090397 

*0082868 

iFUB • 

.5008626 

.1350915 

.0423077 

.0179697 

*0108486 

.0088369 

*0083438 

*0082671 


iTERATISN 17 

9 








J 







K 


40DE 

ri 


12 


13 


14 


15 


16 


17 


18 


QR 



1. 

.23387E 

02 

. 19831E 

02 

. 17563E 

02 

•38142E 

02 

•44243E 

02 

•45356E 

02 

.45543E 

02 

•45562E 

02 

•65906E 

02 

gpi 

\ 1 • 

•23913E 

02 

.20935E 

02 

• 18779E 

02 

•35926E 

02 

*41417E 

02 

.42S08E 

02 

.43301E 

02 

•43498E 

02 

.65928E 

02 

lo< 

21 • 

•23559E 

02 

.21511E 

02 

. 19587E 

02 

.33894E 

02 

.38841E 

02 

•40440E 

02 

•41182E 

02 

•41530E 

02 

•65943E 

02 


31* 

.22723E 

02 

.21649E 

02 

.20032E 

02 

•31991E 

02 

•36456E 

02 

.38209E 

02 

.39158E 

02 

•39634E 

02 

.65955E 

02 

*0t 

41* 

-21610E 

02 

•21429E 

02 

.20159E 

02 

•30178E 

02 

.34220E 

02 

.360S7E 

02 

•37206E, 

02 

.37792E 

02 

.65964E 

02 

■ot 

51* 

•203H4E 

02 

.20916E 

02 

•20010E 

02 

.28428E 

02 

•32099E 

02 

•34050E 

02 

.35310E 

02 

•35991E 

02 

.65972E 

02 

■0( 

61. 

•189H1E 

02 

.20165E 

02 

•19621E 

02 

•26719E 

02 

*30067E 

02 

.32078E 

02 

.33456E 

02 

.34219E 

02 

•65980E 

02 

«0( 

71. 

.17422E 

02 

. 19214E 

02 

. 1902?E 

02 

•25032E 

02 

•28100E 

02 

•30154E 

02 

•31632E 

02 

.32467E 

02 

•659S9E 

02 

■04 

a 1 . 

. 15829E 

02 

. 18090E 

02 

•18232E 

02 

.2334SE 

02 

•26177E 

02 

•28262E 

02 

.29S25E 

02 

.30724E 

02 

•66001E 

02 

|01 

91 • 

.14118E 

02 

.16801E 

02 

•17263E 

02 

•21645E 

02 

•24270E 

02 

•26381E 

02 

.28018E 

02 

.28974E 

02 

.66019E 

02 

To( 

101* 

•12228E 

02 

.15334E 

02 

. 16105E 

02 

•19890E 

02 

.22347E 

02 

*244S3E 

02 

.26189E 

02 

*27igaE 

02 

.66048E 

02 

aQI 

^01. 

•36685E 

00 

• 4600 1 E 

00 

.48316E 

00 

.59669E 

00 

•6704XE 

00 

•73450E 

00 

.78567E 

00 

•81595E 

00 




51* 

•60n7£ 

01 

.29224E 

01 

.22849E 

01 

.22106E 

01 

.20847E 

01 

•19337E 

01 

* 1B261E 

01 

1 17699E 

01 




81 . 

•96851E 

01 

•52S32E 

01 

.41298E 

01 

.38940E 

01 

•35889E 

01 

.32360E 

01 

.29734E 

01 

.28305E 

01 




71* 

• 12339E 

02 

.75111E 

01 

•59691E 

01 

•55985E 

01 

•51372E 

01 

•46005E 

01 

•41914E 

01 

•39647E 

01 




61. 

•14456E 

02 

• 96046E 

01 

•77790E 

01 

.729B1E 

01 

•67035E 

01 

•60027E 

01 

.54579E 

01 

.51517E 

01 




51* 

•16277E 

02 

. 11573E 

02 

.95479E 

01 

.897B2E 

01 

•82716E 

01 

•74265E 

01 

•67576E 

01 

.63770E 

01 




41. 

.17923E 

02 

. 13429E 

02 

. 11270E 

02 

.10630E 

02 

*98306E 

01 

.88599E 

01 

.80789E 

01 

.76294E 

01 




31. 

• 19449E 

02 

. 15182E 

02 

.12939E 

02 

. 12246E 

02 

.11371E 

02 

* 10293E 

02 

.94U8E 

01 

.88992E 

01 




21* 

.20875E 

02 

. 16836E 

02 

.14552E 

02 

.13818E 

02 

•128S4E 

02 

* 11716E 

02 

.10746E 

02 

.10177E 

02 




1 1 * 

•22197E 

02 

. 18390E 

02 

. 16098E 

02 

•15338E 

02 

•14359E 

02 

*13116E 

02 

• 12071E 

02 

• U451E 

02 




1. 

•23387E 

02 

• 19831E 

02 

* 17563E 

02 

.16787E 

02 

•15778E 

02 

•14479E 

02 

•13371E 

02 

•12707E 

02 




^ARAM 

* 1.710344E 

- 1 1 










[ 








huensiqnless 



















^J0DE 

n 


12 


13 


14 


15 


16 


17 


18 


QR 



1. 

. 73473E 

00 

.62300E 

00 

•55174E 

00 

.11983E 

01 

•13899E 

01 

•14249E 

01 

• 14308E 

01 

•14314E 

01 

.65906E 

00 

Sp 

11. 

. 75124E 

00 

.66770E 

no 

•5R995E 

00 

. 11287E 

01 

*13012E 

01 

•13449E 

01 

•13603E 

01 

.13665E 

01 

•65928E 

00 


V 1 • 

• 7 01 4*^ 

no 

, A7i^«of 

no 


00 

• 1 nft48f 

ni 

• t PPOPF 

01 

. 1 ?704P 

Of . 

• 1 

01 

. 1 -■^047F 

Ot 

, AgQ43r 

no 

.n 


31 • 

•71388E 

00 

.68013E 

00 

•62932E 

00 

. 10050E 

01 

•11453E 

01 

. 12004E 

01 

• 12302E 

01 

• 12451E 

01 

• 65955E 

00 *Ak 

41. 

VJ&7S88E 

00 

•67320E 

00 

• 63332E 

00 

•94808E 

00 

•10750E 

01 

. 11337E 

01 

. U68SE 

01 

. 11873E 

01 

•65964E 

00 A 

51. 

.63S50E 

00 

•65710E 

00 

.62864E 

00 

•89310E 

00 

. 10084E 

01 

•10697E 

01 

•11093E 

01 

• U3n7E 

01 

.65972E 

00 

61. 

•59^41E 

00 

•63349E 

00 

•61642E 

00 

•83940E 

00 

•94459E 

00 

•10078E 

01 

•10511E 

01 

•10750E 

01 

.65980E 

00 .01 

71. 

.54733E 

00 

. 60362E 

00 

•59758E 

00 

•78639E 

00 

•88280E 

00 

•94733E 

00 

•99375E 

00 

.lOBOOF 

01 

.65989E 

00 .Bi 

81. 

.49729E 

00 

•56830E 

00 

.57278E 

00 

.73348E 

00 

•82236E 

00 

.887S9E 

00 

.93698E 

00 

•96521E 

00 

•66001E 

00 i|t 

91. 

.44353E 

00 

.52782E 

00 

•54232E 

00 

•67998E 

00 

•76247E 

00 

•82880E 

00 

.88023E 

00 

.91025E 

00 

•66019E 

00 .01 

01. 

.3S417E 

00 

.48172F 

00 

•50596E 

00 

•6248SE 

00 

.70205E 

00 

.76916E 

00 

•S2276E 

00 

.85446F 

00 

•66048E 

00 ^ 

01. 

. 11525E 

•01 

.14451E 

-01 

. 15179E-01 

. 18746E 

-01 

.21062E-01 

.23075E- 

■01 

•24683E- 

01 

•25634E- 

01 



91. 

.18886E 

00 

.91811E 

•01 

.71781E-01 

.69448E 

-01 

•65492E-01 

•60750E-01 

•57370E-01 

•55603E- 

01 



81. 

•30427E 

00 

.16598E 

00 

•12974E 

00 

•12233E 

00 

• U275E 

00 

•10166E 

00 

.93411E- 

01 

. 88923E- 

01 



71. 

•38765E 

00 

•23597E 

00 

•18753E 

00 

•17588E 

00 

•16139E 

00 

•14453E 

00 

•1316BE 

00 

•12455E 

00 



6 1 * 

•45415E 

00 

.30174E 

00 

•24439E 

00 

•22928E 

00 

•21060E 

00 

.18S58E 

00 

.17146E 

00 

•161S4E 

00 



51 . 

•5n37E 

00 

•36359E 

00 

•29996E 

00 

•28206E 

00 

•25986E 

GO 

•23331E 

00 

•21230E 

00 

•20034E 

00 



41. 

.5630SE 

00 

•42189E 

no 

•35404E 

00 

•33395E 

00 

.308S4E 

00 

•27834E 

00 

•25381E 

00 

•23968E 

00 



31. 

•61 102E 

00 

•47695E 

00 

•40650E 

00 

.38471E 

00 

•35724E 

00 

.32337E 

00 

•29568E 

00 

-27958E 

DO 



21. 

.655S2E 

00 

.52S92E 

00 

•45715E 

00 

.43412E 

00 

•40477E 

00 

•36806E 

00 

•33761E 

00 

.31971E 

00 



11. 

•69735E 

00 

•57773E 

00 

•50573E 

00 

•48185E 

00 

•45109E 

00 

•41207E 

00 

•37922E 

00 

•35973E 

00 



1. 

•73473E 

00 

.62300E 

00 

•55174E 

00 

.52739E 

00 

•49569E 

00 

.454SSE 

00 

.42006E 

00 

•39919F 

00 



EFLECTANCE FIELD 

















TERATI0N 4 








0DE 11 

12 

13 

14 

15 

16 

17 

18 


1. 

•31385E 

02 

.41866E 

02 

•44539E 

02 

.45277E 

02 

•45490E 

02 

•45549E 

02 

•45561E 

02 

•45562E 

02 

11. 

•30230E 

02 

•39166E 

02 

•41734E 

02 

.42672E 

02 

.43104E 

02 

•43334E 

02 

.43462E 

02 

•43524E 

02 

21. 

.2908SE 

02 

•36717E 

02 

•39170E 

02 

•40259E 

02 

•40866E 

02 

•41238E 

02 

•4146AE 

02 

•415S0E 

02 

31. 

•27931E 

02 

.34456E 

02 

.36792E 

02 

•37993E 

02 

•38742E 

02 

•39234E 

02 

•39542E 

02 

•39704E 

02 

41. 

• 26744E 

02 

.32338E 

02 

.3455SE 

02 

•35S42E 

02 

.36707E 

02 

.37298E 

02 

.37678E 

02 

.37879E 

02 

51. 

■25517E 

02 

•30326E 

02 

.32436E 

02 

•33782E 

02 

•34741E 

02 

•35416E 

02 

•35857E 

02 

•36093E 

02 

61. 

•24246E 

02 

•28394E 

02 

.30402E 

02 

•31792E 

02 

•32828E 

02 

•33575E 

02 

.34069E 

02 

•34335E 

02 

71. 

•22925E 

02 

•26518E 

02 

•28432E 

02 

•29854E 

02 

•30953E 

02 

•31761E 

02 

.32302E 

02 

.32595E 

02 

81. 

•21546E 

02 

•24674E 

02 

•26504E 

02 

•27950E 

02 

.29102E 

02 

.29963E 

02 

.30545E 

02 

•30862E 

02 

91. 

•20098E 

02 

.22836E 

02 

.24594E 

02 

•26059E 

02 

•27257E 

02 

•28165E 

02 

•28784E 

02 

•29122E 

02 

01. 

• 18554E 

02 

•20970E 

02 

•22668E 

02 

•24151E 

02 

•25391E 

02 

•26343E 

02 

•26996E 

02 

•27355E 

02 


01* 

* 5566HE 

00 

•62909E 

00 

•68004E 

00 

•72454E 

91.* 

.•22181E 

01 

•21711E 

01 

•20619E 

01 

•19564E 

8l'. 

•39430E 

01 

•37910E 

01 

•35362E 

01 

•32899E 

71. 

•56826E 

01 

.54408E 

01 

•50578E 

01 

.46833E 

61. 

•74113E 

01 

.70944E 

01 

•66008E 

01 

•61119E 

51. 

.91154E 

01 

•87366E 

01 

.S1488E 

01 

•75593E 

41, 
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•71283E»0l 

*71283E-01 

.71283E-01 

•71283E-01 

.712S3E-01 

• 71283E-.01 

•71283E-01 

.71283E-01 

•71283E-01 

•712S3E 

•712S3E-01 

•712S3E-01 

.712S3E-01 

•71283E-01 

*71283£**01 

•71283E-01 

.71283E-01 

•71283E-01 

.71283E-01 

• 712 fl 

E 

•712S3E-01 

.712S3E-01 

•712S3E-01 

•71283E-01 

*71283E-01 

*71233E-01 

-71233E-01 

.712S3E-01 

.71283E-01 

•71283E 

.71283E-01 

• 712R3E-01 

.71233E-01 

•71283E-01 

•71283E-01 

.71283E-01 

.71283E-01 

.712832-01 

•71283E-01 

* 712 ! 

E 

.71283E-01 

•71283E-01 

.71283E-01 

•71283E-01 

•71283E-01 

•71283E-01 

*71283E-0l 

•71283E-01 

.71283E-01 

*7X2M 

E 

•712S3E-01 

•71283E-01 

•71283E-01 

. 71283E-01 

•712S3E-01 

.71283E-01 

•71283E-01 

•712S3E-01 

.71283E-01 

* 712^ 

E 

• 71283E-01 

•71283E-01 

•712S3E-01 

.712S3E-01 

•71283E-01 

• 71283.E-01 

•71283E-01 

.71283E-01 

.71283E-01 

.712S3E-01 

*712^ 

E 

•71283E-01 

•71283E-01 

*71283E-01 

•712S3E-01 

.71283E-01 

.71283E-01 

*71283E-01 

i ’ >■ ^ T 

• 71283E-01. 

•71283E 

•71283E-01 

•71283E"01 

•7l?83E-01 

*7I2S3E-01 

•712a3E-01 

-712S3E-01 

-71283E-01 

*71283£-01 

.71283E-01 

*7l2j 

E 

•71283E-01 

"71283E-01 

. 71283E-01 

•712S3E-01 

•71283E-01 

•71283E-01 

*71283E-01 

.71283E-01 

•71283E-01 

*712£tf 



•71283E-01 
CRT = 1*57080 
CRTOG “ S9*y99y 
RFLI «■ *0000000 

*0000000 

*0000000 

* ooooooo 

•ooooooo 

.ooooooo 

1 

•ooooooo 

•ooooooo 

RFL0 * *0000000 

•OOOOOCO 

* ooooooo 

•ooooooo 

•ooooooo 

.ooooooo 

•OOOOOOO 

•ooooooo 


ITERATION ■ 28 


- 1 1 ■' 

12 


13 


14 


15 


16 


17 


T8 


QR 


6qro| 


10429E 

02^^.t0429E 

0.2 

• 10429E 

02 

.10429E 

02 

. 10429E 

02 

. 10429E 

02 

. 10429E 

02 

. 10429E 

02 

•45166E 

01 

-•19543E 

C 

94S92E 

01 / .98125E 

01 

• 99S44E 

01 

• 10076E 

02 

.10130E 

02 

•10362E 

02 

. 10182E 

02 

• 10191E 

02 

. 44842E 

01 

-.18671^ 

C 

•91077E 

01 .93593E 

01 

.957032 

01 

.97036E 

0: 

.9799SE 

01 

•9S5S0E 

01 

.9S937E 

01 

•99121E 

01 

•44630E 

01 

-il7948| 

c 

• 8755f)E 

0l'^.39754E 

01 

.91S69E 

01 

. 93479E 

01 

•94613E 

01 

.95380E 

01 

.9586^E 

01 

•961 Igr 

01 

• 44440E 

01 

-•17252f 

c 

.84107E 

OUU. 86183F. 

01 

•S8227E 

01 

.S9927E 

01 

.91202E 

01 

.92099E 

01 

•92679E 

01 

.92987F 

01 

•44262E 

01 

-•16566En?o 

•80692E 

01§.«271SE 

01 

•R4699E 

01 

•86430E 

01 

.S7788E 

01 

.88774E 

01 

.89426E 

01 

•89776R 

01 

•44094E 

01 

••158S5| 

>0 

•77236E 

01^* 79294E 

01 

.S1234E 

01 

•82969E 

01 

.S4376E 

01 

•85424E 

01 

•86128E 

01 

•865tie 

01 

•43937E 

01 

-•152051 

^0 

• 73B6>6E 

01 .75875E 

01 

.77794E 

01 

•79528E 

01 

.80963E 

01 

.82054E 

01 

.S2797E 

01 

*832032 

01 

•43791E 

01 

14520E 

0 

• 70390E 

013. 72425E 

01 

.74343E 

01 

.76078E 

01 

.77534E 

01 

•7S656E 

01 

•7942SE 

01 

.79S54E 

01 

•43661E 

01 

-•13821ft 

0 

66766E 

01^.68876E 

01 

. 70S24E 

01 

.72575E 

01 

•74052E 

01 

.75201E 

01 

.75998E 

01 

•764412 

01 

•43561E 

01 

-•130821 

0 

■-62416E 

OlJ *649072 

01 

.67002E 

01 

•68830E 

01 

.70360E 

01 

•71552E 

01 

•72382E 

01 

•72843E 

01 

•43662E 

01 

-.12078E 

0 


-■49933E 

OU • 

51926E 

01 

•53602E 

01 

•55064E 

01 

•56288E 

01 

•57241E 

01 

•57905E 

01 

•58275E 

01 

-64198E 

Oil. 

60091E 

01 

•5S646E 

01 

.58468E 

01 

•58745E 

01 

•59131E 

01 

•59465E 

01 

•59669E 

Cl 

63246E 

01 1 . 

65389E 

01 

•63207E 

01 

•62166E 

01 

•61788E 

01 

•61723E 

01 

.61772E 

01 

•61B28E 

01 

- 71797E 

01^. 

69485E 

01 

•67267E 

01 

•65819E 

01 

•65020E 

01 

• 64626E 

01 

•64451E 

01 

• 64386F 

01 

•75244E 

01^* 

7313SE 

01 

•71019E 

01 

•69392E 

01 

.68333E 

01 

•67703E 

01 

•67357E 

01 

•67196E 

01 

-78655E 

Olu.* 

76620E 

01 

•74598E 

01 

•72901E 

01 

•71682E 

01 

•70886E 

01 

•70411E 

01 

•70174E 

01 

•S2062E 

oif. 

80050E 

01 

•780S9E 

01 

.76363E 

01 

•75050E 

01 

•74138E 

01 

•73566E 

01 

-73272E 

01 

•S5487E 

01'^. 

83474E 

01 

•81543E 

01 

•79815E 

01 

.78437E 

01 

•77442E 

01 

•76798E 

01 

•76459E 

01 

-■38957E 

Olii. 

86924E 

01 

• 85000E 

01 

•83267E 

01 

•81848E 

01 

•80795E 

01 

•80096E 

01 

•79721E 

01 

•92527E 

OIT* 

90443E 

01 

•88502E 

01 

•86757E 

01 

•85305E 

01 

•84206E 

01 

.83464E 

01 

•83061E 

01 

•96517E 

Old • 

94213E 

01 

•92185E 

01 

•90392E 

01 

•88896E 

01 

•S7750E 

01 

•86967E 

01 

•86539E 

01 


» 6.465228E-11 
p/h 


'0NLESS 

S>f/^CAJs/i:>A^i.Bss 


11 

12 


13 


14 


15 


16 


17 


18 


OR 


DORojl 


ICOOOE 

01^ . lOOOOE 

01 

• lOOOOE 

01 

• lOOOOE 

01 

. lOOOOE 

01 

•lOOOOE 

01 

•lOOOOE 

01 

• lOOOOE 

Cl 

.13785E 

00 

-.37562E 

0 

90938E 

00] .94087E 

00 

•95736E 

00 

•96615E 

00 

•97129E 

00 

•97442E 

00 

.97627E 

00 

•97721E 

00 

. 13686E 

00 

-•36886|| 

0 

87330E 

Ooj .89742E 

DO 

•91766E 

00 

•93101E 

00 

•93966E 

00 

•94524E 

00 

.94S66E 

00 

'•95043E 

00 

•13622E 

00 

-•344971 

0 

83952E 

0C«4*86061E 

00 

•8S089E 

00 

•89633E 

00 

•90720E 

00 

•91455E 

00 

•91920E 

00 

•92163E 

00 

• 13564E 

00 

-•33157E 

0 

80646E 

00 .82637E 

00 

•84597E 

00 

•86227E 

00 

.S7449E 

00 

•S8309E 

00 

tS8S66E 

00 

•89161E 

00 

•13509E 

00 

-•31S40g 

^0 

77372E 

OoS*79314E 

00 

•81214E 

00 

•82873E 

00 

•84i75E 

00 

•85121E 

00 

.85746E 

00 

• 860SaE 

no 

• 13458E 

00 

-•3053li 

»0 

74106E 

00®. 76032E 

00 

•77892E 

00 

•79556E 

00 

•80904E 

00 

•81909E 

00 

.82584L 

00 

•82951E 

00 

•13410E 

00 

-.29224?T 

^0 

7CS26E 

00 .72754E 

00 

•74593E 

00 

•76255E 

00 

•77632E 

00 

•78678E 

00 

•79390E 

00 

•79780E 

no 

•13365E 

00 

-.27908^ 

0 

67494E 

00^.69445E 

00 

•712S5E 

00 

•72948E 

00 

•74344E 

00 

.75419E 

00 

•76160E 

00 

•76568E 

00 

•13326E 

00 

-•265641 

0 

64019E 

00>'.66042E 
00<| .62237E 

00 

•67910E 

00 

•69588E 

00 

•71005E 

00 

•72107E 

00 

•72871E 

00 

•73295E 

00 

•13295E 

00 

-•25143P 

0 

59848E 

00 

•64245E 

00 

•65998E 

00 

•67465E 

00 

•6860SE 

00 

•69403E 

00 

•69S46E 

00 

•13326E 

00 

-•23213E 

0 


01 . 

•47S79E 

00 

.49789E 

00 

•51396E 

00 

.52798E 

00 

91* 

•61S56E 

00 

.57619E 

00 

•56233E 

00 

•56062E 

00 

81,. 

•65438E 

00 

.62699E 

00 

.60607E 

00 

.59608E 

00 

71. 

.68843E 

00 

. 66626E 

00 

.64499E 

00 

•63110E 

00 

61* 

.72148E 

00 

.70125F 

00 

•68097E 

00 

t66537E 

00 

51. 

•75419E 

00 

•73468E 

00 

•71529E 

00 

•69902E 

00 

41. 

•78685E 

00 

•76757E 

00 

.74876E 

00 

.73226E 

00 

31. 

•81969E 

00 

.80039E 

00 

*7S18SE 

00 

•76531E 

00 

21 . 

.85297E 

00 

.83347E 

00 

.81503E 

00 

.79841E 

00 

u,. 

•38720E 

00 

.86722E 

00 

•84861E 

00 

•831S7E 

00 

1. 

•92546E 

00 

.90336E 

00 

•88392E 

00 

•86672E 

00 


R(KK) « .862564 


•53972E 00 .54S86E 00 
•56328E 00 •S6698E 00 
• 5924SE 00 .?i9l84E 00 
•62345E 00 •61967E 00 
•65521E 00 •64917E 00 
•68732E 00 .67969E 00 
t7l962E 00 .710S7E 00 
•75210E 00 .74256E 00 
•78480E 00 .77470E 00 
•81795E 00 .80741E 00 
•85238E 00 .84140E 00 


. 55523 E 00 .558/7E 00 
.5701SE 00 .572t3r 00 
•59231E 00 .59a84E 00 
•61799F 00 -61737E 00 
•64586E 00 .64431E 00 
•67514E 00 *67287^ 00 
•70539E 00 .70257E 00 
.73639E 00 .73313E 00 
•76800E 00 .76441E 00 
.80030E 00 .79643E 00 
.83389E 00 .82978E 00 


« O t_ -j ’ 


0VERALL FLUXES> STARTING AT FR0NT WALUS 


-63371E 

02 

•6329SE 

02 

.63237e 

02 

•63184E 

02 

•62916E 

02 

.628S4E 

02 

•62S52E 

02 

.62822E 

02 

•62619E 

02 

.62591E 

02 

• 62563E 

02 

• 62536E 

02 

.62351E 

02 

f 6232SE 

02 

.62299E 

02 

.62274E 

02 

•62100E 

02 

• 62076E 

02 

•62052E 

02 

•62028E 

02 

•61864E 

02 

.61841E 

02 

.61818E 

02 

.61795E 

02 

.61641E 

02 

•61620E 

02 

.61599E 

02 

.61578E 

02 

.61435E 

02 

•61416E 

02 

•61396E 

02 

.61377E 

02 

•61HblE 

02 

•61235E 

02 

.61219E 

02 

.612Q3E 

02 

•61111E 

02 

.61 102E 

02 

.61095E 

02 

• 61090E 

02 


•6ia51£ 02 

0VERALL divergence: 

-•27A17E 01 -•27260E 01 -.27117E 01 -i26985E 01 

-•26193E 01 “'adOSOE 01 -«*?5985E 01 -•25S83E 01 

-•25179E 01 ••250S0E 01 -.249S2E 01 -t24883E 01 

-*24202E 01 -•2A105E 01 “«24008E 01 -t239l2E 01 

-»23240£ 01 "•23144E 01 -•23048E 01 -•22953E 01 

-•222S5E 01 -<22189E 01 -*22094E 01 -*21999E 01 

-•21330E 01 01 "•2U39E 01 -•21043E 01 

-•20370E 01 -•20273E 01 “i20l76£ 01 -*20078£ 01 

-•193S9E 01 -*19238E 01 -•19188E 01 -*190S6E 01 

-•18351E 01 -•18240E 01 -•18127E 01 -»1S010£ 01 


•63137E 

02 

• 63094E 

02 

fl 63055E 

02 

.63018E 

02 

• 62982E 

.62792E 

02 

• 62762E 

02 

•62733E 

02 

.62704E 

02 

.62675E 

•62509E 

02 

• 62482E 

02 

. 62456E 

02 

•62429E 

02 

•62403E 

.62249E 

02 

• 62223E 

02 

.62198E 

02 

.62174E 

02 

.62149E 

.62004E 

02 

• 61980E 

02 

•61957E 

02 

• 61933E 

02 

• 61910E 

.6X773E 

02 

• 61751E 

02 

.61729E, 

02 

.61707E 

02 

•61685E 

•61557E 

02 

.61536E 

02 

•61516E 

02 

• 6149bE 

02 

•61475E 

•61359E 

02 

•61340E 

02 

• 61322E 

02 

i6l304E 

02 

.61286E 

#611S7E 

02 

•61173E 

02 

.61159E 

02 

.61145E 

02 

• 6U33E 

.61089E 

02 

-61092E 

02 

.61100E 

02 

•611 17E 

02 

.61144E 


.26S60E 01 -.H67^2E 01 -.26627E 01 -.2691feE 01 -.26406E^1 
'.25781E 01 -.25679E 01 -.2557SE 01 ".2547SE 01 -.25378 e1i 
. 24785E 01 -.24688E 01 -.24590E 01 -.24493E 01 -.24395EJ1 
'.23816E 01 -•23719E 01 -.23623E 01 ••23431E 01 
-.22857E 01 -.22762E 01 -.22666E 01 -.92^7,1E 01 -*22475Ejl 
'•21903E 01 -.21808E 01 -.21712E 01 -.21617E 01 -.21522EJ1 
..20947E 01 -.20851E 01 --20755E 01 -.?0659E 01 -.20563E^: 
■.19381E 01 -.19883E 01 -.19785E 01 -.19686E 01 -.195S7EJ1 
».18984E 01 -fl8882E 01 -.18778E 01 -.1R673E 01 -.18567E 01 
•.17839E 01 -*17762E 01 -.17628E 01 -.17483E 01 -.17325Eft 



Tes'r Oa.se, 3 

^^=•5.(77 


0 "K 

X^O S«^^Tr£/i//u<il 

Jj01>£X - r C ay 

Sl/\B Ttoo LAV£.RS — ou~ /S’ 


NDS - 201 

T8L » 9.999999E-11 
TEST - 1000 
PRINT » 25 
N0NDM B 1 

THICK =' 1.00000 CH 
BNDS - 2000000000 

SPECIAL — 1 STANDARD 
TEMPO *1 

TEMP » .000 .000 .000. . .000 *000 

IKELVINI 


to " *^99^ Re/\>a. ^ 


.000 .000 *000 tOOO .000 



isar = 


• 000 

• 000 

• 000 

• 000 

• 000 

• 000 

• 000 

• 000 

•000 

• 000 

• 000 

•000 

• 000 

• 000 

• 000 

*000 

• 000 

• 000 

• 000 

.000 

• 000 

• 000 

•000 

• 000 

• 000 

• 000 

• 000 

..000 

• 000 

.000 

• 000 

• 000 

• 000 

• 000 

• 000 

.000 

.000 

• 000 

• 000 

.000 

• 000 

• 000 

• 000 

• 000 

• 000 

• 000 

.000 

.000 

• 000 

.000 

• 000 

• 000 

• 000 

• 000 

• 000 

.000 

.000 

• 000 

.000 

•000 

• 000 

• 000 

•000 

• 000 

• 000 

.000 

• 000 

• 000 

•000 

• 000 

• 000 

• 000 

• 000 

• 000 

• 000 

.000 

.000 

• 000 i 

• 000 

.000 

• 000 

• 000 

• 000 

• 000 

• 000 

.000 

• 000 

• 000 

• 000 

.000 

• 000 

• 000 

♦ 000 

• 000 

.000 

.000 

.000 

.000 

• 000 

.000 

• 000 

• 000 

• 000 

• 000 

• 000 

.000 

.000 

• 000 

• 000 

.000 

.000 

• 000 

• 000 

• 000 

.000 

• 000 

.000 

• 000 

• 000 

.000 

.000 

• 000 

• 000 

• 000 

.000 

.000 

.000 

• 000 

• 000 

.000 

• 000 

• 000 

• 000 

• 000 

.000 

.000 

.000 

.000 

• 000 

• 000 

.000 

• 000 

• 000 

• 000 

• 000 

• 000 

.000 

• 000 

• 000 

• 000 

• 000 

• 000 

•000 

• 000 

.000 

.000 

.000 

• 000 

•000 

.000 

• 000 

• 000 

• 000 

• 000 

.000 

• 000 

.000 

.000 

• 000 

.000 

• 000 

• 000 

• 000 

• 000 

.000 

.000 

.000 

• 000 

.000 

.000 

.000 

0000000 

•0000000 

•0000000 • 

0000000 

•ooooooo 






i 


INTENSITIES IN W ATTS/CM*^2/STERA0I AN' 
FLUXES IN WATTS/CM4»^2 


FLUX DIVERGENCE IN WATTS/CM¥¥3 



^ in 


SPECIAL BAND 


BAND 1 


INDEX « 1*000 

RB « *800 

Q0 » 100*0000000 WATTS/CM^f^a 

» *158850 *158850 *001588 *001588 *001588 CM¥#-1 

^ 3*018150 3*018150 3*175412 3*175412 3*175415 


GAUSSIAN integral 

MU'S WEIGHTS 

•09501249 # 18945062 

•28160357 *18260342 

• 45801675 * 16915649 

•61787623 *14959598 

• 75540441 • 12462896 

*86563122 *09515852 

•94457501 *06225352 

•98940092 *02715246 


WAVELENGTH INTERVAL; *OOOOOE 00 - *99999E 70 CM 

•OOOOOE 00 • *99999E 74 MICRONS 


F(0 - LT) : 


• 76529E- 

•80 

•41180E-82 .10295E-S3 

•51476E-84 

•OOOOOE 00 

•OOOOOE 

00 

•OOOOOE 00 .OOOOOE 00 

•OOOOOE 00 

•OOOOOE 00 

• 15674E 

32 

•25098E 00 *OODOOE 00 

•OOOOOE 00 

•OOOOOE 00 

• 32680E' 

•79 

•26313E-77 *62500E-01 

•OOOOOE 00- 

•65446E-77 

•OOOOOE 

00 

•OOOOOE 00 *43959E"72 

•86736E-18 

•22041E-38 

•22S59E 

09 

•40116E-17 .27636E-75 

•27537E 48 

•44060E 49 

•27538E 

48 

•23510E-37-.96207E 12 

•78603E-21 

•27638E-75 

•35112E 

51 

•43656E-10 *11755E-37 

•36734E-39 

•l2l4aE-26 

•33144E 

13 

•35903E-73 .36734E-39 

.18808E-36- 

.21690E-06 

*6o8i8E 

08 

•25760E 51 *2754lE 48 

•94039E-37 

•20881E-52 

,r?97i Hr 

40 

• nntinn'“ nn F ®n 

• 1 1 1 9'^'“ 



• OOOOOE 

00 

•52964E-77 

•21620E-81 tOOOOOE 00 

•51476E 

* OOOOOE 

00 

•OOOOOE 00 

•OOOOOE 00 •OOOOOE 00 

• 000 


•OOOOOE 

00 

•OOOOOE 00 

•OOOOOE 00 •51476E-S4 

• 658 

|,E 

•OOOOOE 

00 

•80964E-77 

•22204E-15 *13092E-79 

•503 

T5E 

• 46137E 

00 

0 

• 35112E ^51 

•13605E 46 •25759E 51 

•188 


• 94o4oE* 
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